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Abstract Over the last two decades the uninterrupted, high resolution observa- 
tions of the Sun, from the excellent range of telescopes aboard many spacecraft 
complemented with observations from sophisticated ground-based telescopes have 
opened up a new world producing significantly more complete information on the 
physical conditions of the solar atmosphere than before. The interface between the 
lower solar atmosphere where energy is generated by subsurface convection and 
the corona comprises the chromosphere, which is dominated by jet-like, dynamic 
structures, called mottles when found in quiet regions, fibrils when found in active 
regions and spicules when observed at the solar limb. Recently, space observations 
with Hinode have led to the suggestion that there should exist two different types 
of spicules called Type I and Type II which have different properties. Ground- 
based observations in the Can H and K filtergrams reveal the existence of long, 
thin emission features called straws in observations close to the limb, and a class 
of short-lived events called rapid blue-shifted excursions characterized by large 
Doppler shifts that appear only in the blue wing of the Can infrared line. It has 
been suggested that the key to understanding how the solar plasma is accelerated 
and heated may well be found in the studies of these jet-like, dynamic events. 
However, while these structures are observed and studied for more than 130 years 
in the visible, but also in the UV and EUV emission lines and continua, there 
are still many questions to be answered. Thus, despite their importance and a 
multitude of observations performed and theoretical models proposed, questions 
regarding their origin, how they are formed, their physical parameters, their asso- 
ciation with the underlying photospheric magnetic field, how they appear in the 
different spectral lines, and the interrelationship between structures observed in 
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quiet and active regions on the disk and at the hmb, as well as their role in global 
processes has not yet received definitive answers. In addition, how they afltect the 
coronal heating and solar wind need to be further explored. In this review we 
present observations and physical properties of small-scale jet-like chromospheric 
events observed in active and quiet regions, on the disk and at the limb and discuss 
their interrelationship. 

Keywords Sun ■ chromosphere • spicules ■ jet-like structures ■ small-scale events • 
physical parameters 



1 Introduction 

During the last 2 — 3 decades, observations from space combined with ground- 
based data, as well as state-of-the-art, post-processing reconstruction techniques 
combined with progress in numerical modelling and theory have increased our 
understanding of the physical processes occurring in the solar atmosphere. Apart 
from demonstrating that the solar atmosphere is highly dynamic and inhomoge- 
neous these studies have also provided a deep insight into all phenomena occurring 
in its different parts. They have revealed a number of different types of structures 
that exist in the different atmospheric layers for which a range of formation mech- 
anisms have to be investigated. They have also provided important information 
about the physical conditions of the emitting plasma, as well as hints about the 
physical nature of related dynamic phenomena which continuously occur on the 
Sun over a large range of spatial and temporal scales. 

The solar chromosphere is an intricately structured and dynamic layer of the 
Sun. It is traditionally defined as the layer, about 2 000 km thick, lying between 
the photosphere and the transition region and has long been known not to be in 
hydrostatic equilibrium. Therefore, its nature is that of dynamic fine-scale struc- 
tures which uniquely characterize this part of the solar atmosphere. These appear 
to be structured by the magnetic fields that extend from network boundaries and 
plages and be directly related to its evolution. The appearance of the underlying 
atmosphere is, mainly, due to processes involving the emergence, buffeting and 
re-arrangement of the magnetic fiux generated by interior dynamo action. The 
resulting spatial inhomogeneity produces a fruitful environment for the emerging 
fiuxes to undergo quite complex changes as they dynamically evolve and inter- 
act with the pre-existing magnetic fields in the ambient medium. Complex and 
concentrated, mainly unipolar, strong magnetic fields produce the so-called active 
regions, which are the locations of some jet-like structures, such as the (traditional) 
fibrils and dynamic fibrils (DFs). Outside active regions, what is generally known 
as the quiet solar chromosphere, is not really quiet. It is continually subject to a 
variety of small-scale events, such as mottles, straws and the, recently reported, 
rapid blue excursions, occurring at the boundaries of cellular patterns that con- 
stitute the magnetic network. The magnetic network is believed to be sustained 
by magnetic bipoles, which emerge in the interiors of supergranules, move apart 
and are fragmented by granular buffeting. Eventually, the fragments are driven 
to t he supergarnular boundaries where they interact with the existing magnetic 
fiux (jWang et al.l Il996t ISchriiver et al.l Il997l) . Opposite polarity fiuxes cancel and 
submerge, whereas like-polarity flux merges to form larger flux concentrations. 
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When observing the chromosphere at different wavelengths, from the EUV 
to the radio domain its very large range of inhomogeneity and structuring is re- 
vealed. Chromospheric structures can be seen against the solar disk, by means of 
monochromatic or narrow-band filters or spectrometers operating in strong spec- 
tral lines, such as Ha, Can H & K, and the Can infrared triplet (IR). In particular, 
filtergrams recorded in the wings and line center of Ha reveal a wealth of fine-scale 
structures. The most prominent small-scale features residing at the network bound- 
arie s are ce r tainly mottles. Mottles are thin jet-like features, better observed in 
Hg. lRuttenl ( 2007|) identified in Can H images taken close to the solar limb as very 
thin and short-lived bright structures which occu r in "he dge rows " . He called them 
"straws" . Recently, Langangen et al. (2008^ and lRoupp e van der Voort et al. 
found sudden, large line shift changes in both Ha and Can IR lines on their blue 
side. The authors interpreted these events, which occur at the edges of rosettes 
(cluster of mottles expanding radially around a common center) , as chromospheric 
up- flows and called them "rapid blue-shifted excursions" (RBEs). In active re- 
gions, a subset of dark structures called fibrils includes relatively short, thin, jet- 
like structures called dynamic fibrils (DFs). They are found in the central area 
of plages where the magnetic field is more dense, and more vertically oriented. 
Another large subset includes very inclined, almost horizontal, long fibrils which 
are clearly more stable than the DFs. The relationship between all these on-disk 
structures, as well as their relationship to spicules, thin, elongated, jet-like struc- 
tures observed at the solar limb, is under strong debate. Spicules are the most 
obvious component of the chromospheric limb. Although they were discovered al- 
most 130 years ago, they still remain one of the mysterious phenomena in the solar 
atmosphere because, at least until recently, remained very close to the resolution 
limits of current solar observations. As a result, several of their properties have 
not been well defined and this h as led to a pleth ora of theoretical models for their 
interpretation (for a review see ISterlind ( 20001) 1. Recently, with its high resolu- 
tion and high cadence, as well as an uninterrup ted view of the Su n, the Hinode 
space mission provided a new picture of spicules. iDe Pontieu et ahl ([2007b) , based 
on time series observations of the solar limb taken with Solar Optical Telescope 
(SOT) aboard Hinode in the Can H broad filter, suggested that spicules can be 
grouped in two categories, called Type I and Type II, which seem to have dif- 
ferent properties and formation mechanisms. The origin of short duration/high 
velocity spicules (Type II) has been attributed to magnetic field reconnection. On 
the other hand, the origin of Type I, classic spicules (long duration/low velocities) 
seems to be associated with the leakage of p-modes into the upper atmosphere and 
subsequent development of shocks that follow the magnetic field lines. 

While spicules are ubiquitous at the li mb, there is still a debate on the i s sue o f 
their counterparts are on the solar disk. iGrossmann-Doerth and Schmidt! (|1992|) 
concluded from the different velocity distributions of m ottles and spicules, that 
mottl es are not the disk counterparts of (classical) spicules. iTsiropoula and Schmieden 

on the other hand, based on the similarity of their physical parameters, ar- 
gued for a close relationship between them. Type I spicules appear to rise up from 
the limb a nd fall back again. They s how a similar dynamical behaviour as active re- 
gion DFs ( De Pontieu et aLLl2007a|) . as well as a subset of mottles. Type II spicules, 
on the other hand, seem to exhibit an upward motion followed by rapid fading from 



the Hinode Ca ii H passband, although more recent work bv [ Zhang et al 



201S 



were unable to locate any Type II spicules. iRouppe van der Voort et al.l (j200 
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suggested that the disk counterpart of Type II spicules are RBEs based on the 
similarities of their properties. 

When observed in the line center of Ha, the chromospheric plasma is clearly 
seen to be organized along fine-scale magnetic stru ctures which have different 
inclinations. Recentlv. lKontogiannis et al.l ( 2010al) andlKontogiannis et al.l ( 2010b|) 



using high-resolution Ha observations together with magnetograms obtained with 
the spectropolarimeter on Hinode, revealed through a potential magnetic field 
extrapolation another very important aspect of these structures when observed on 
the disk. They showed that magnetic flux tubes following the local inclination of 
the magnetic field lines, define the layer of the magnetic canopy, i.e. the plasma- 
P I layer, where the plasma-/? is the ratio of gas pressure to magnetic pressure. It 
is at this layer, where acoustic waves generated by the convective motions, undergo 
mode conversion, refraction, reflection and transmission with important effects to 
the upper solar atmosphere and thus these structures play a very important role 
in wave propagation. 

The fact that the solar corona has a temperature more than 1 million degrees 
is still a puzzling problem of solar physics, despite the considerable theoretical 
and experimental efforts. Since the energy release in the largest explosive events 
(flares and microflares) does not supply enough power to heat the corona, the 
behaviour of smaller-scale and less e nergetic, but much more frequent event s is an 
essential constituent of the problem ( Dovle and Butleil . ll985HParkeiill988l) . Both 



magnetic reconnection and waves related to spicules are powerful mechanisms for 
energy release and transfer and have been considered as hea ting agents of the solar 
corona. pTsiropoula et al.l ( 1994h and iTziotziou et~al ] (|2003l) examining the tempo- 



ral evolution of the line-of-sight velocity of mottles found bi-directional flows and 
proposed a model according to which spicules are due to magnetic reconnection. 
iTsiropoula and Tziotzioul ( 20041) showed that if magnetic reconnection is consid- 



ered as the driving mechanism of mottles, the material they supply to the solar 
corona is in excess of that needed to compensate for coronal mass losses in the 
solar wind, while the amount of energy released to heat the corona depends on 
several parameters, among which are the number of the events, t heir axial velocity 
and m agnetic field, etc., and can be negligible or substantial. iDe Pontieu et al.l 
on the other hand, showed that the energy flux provided by the Alfven 
waves identified in spicules is large enough to supply the energy necessary to 
heat the quiet solar corona and drive the solar wind. Thus the question arises 
of how these small-scale features participate and contribute into the larger scale 
observable phenomena, and how a more accurate determination of the physical 
parameters, which are related to the physical processes, can be extracted from 
observations. 

This review is the first part of a series of four reviews. The aim of these reviews 
is to present an overview of our current understanding (observations, modelling, 
physical parameters, etc.) of solar small-scale, jet-like structures observed both 
on the disk and at the limb, in active, as well as in quiet regions. In this review 
(Part I) we present an overview of the chromospheric fine-scale structures. More 
specifically, we will present their morphological properties, their derived physical 
parameters, such as height, width, lifetime, velocities, etc., as well as methods 
used for the derivation of some of them, their dynamical behaviour, and comment 
on their interrelationship. In subsequent papers we will discuss transition region 
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small-scale events (Part II) , and the role of the various small-scale events in coronal 
heating (Part III). 



2 The solar photosphere 

The solar surface is threaded by a complex network of magnetic fields. Observa- 
tions of these fields at the photospheric level reveal them to be concentrated into 
regions of mixed (negative and positive) polarities covering a wide range of values 
for the (unsigned) magnetic flux and lifetimes. In quiet areas the dynamic inter- 
play between the magnetic fie ld concentrations and the plasma has been called 
the Sun's "magnetic carpet" ( Title and Schriiveil Il998h . In these areas, photo- 



spheric motions are dominated by the fiow pattern of large convective cells called 
supergranules. These cells range in diameter from roughly 10 OOP km to 50 O OP km , 
with an average diameter of between 13 000 and 18 000 km (|Hagenaar et aHligQ?!) . 
Their flow pattern takes the form of an up-flow at the cell cen ter, followed by radial 
outflow of ~0.5kms~^, and do wn-flow at the cell boundaries ( Simon and Leightonl 



Il964|) . IWang and ZirirJ (11989 ) derived an upper limit of O.lkms"^ for the mag- 
nitude of the supergranule's vertical velocity (both up-flow and down-flow). The 
strongest down-flow occurring at points where two or more cells meet, and as a 
result, the magnetic flux tends to build up along the boundaries of supergranu- 
lar cells after being swept from the center by the radial outflow. It is now well 
established from investigations of high resolution magnetograms that new bipo- 
lar elements emerge continuously inside the cell interiors and are, subsequently, 
swept by the supergranular flow at the periphery of supergranular cell s where they 
gather and form the netwo rk boundaries with kG field strengths (IWang et al.l 
Il996l: ISchriiver et~aD . Il997h . This is illustrated weU by the SoHO/MDI magne- 
togram shown in Figure[T] The supergranule cell boundaries are marked in yellow, 
and it can be seen that the majority of the flux concentrations are located along 
these lines, and in particular at the intersection of multiple cells. These concentra- 
tions called network concentrations are typi cally found at sites of strong down-flow, 
at the boundaries of supergranular cells (Martin, '1988). Their fluxes ar e of the 
order of 10^** - 10^^ Mx and have typical diameters of 1 000 - 10 000 km (iParnelll 
I2OOII) . They tend to be larger than ephemeral regions which are clumps of newly 
emerged flux of varying strengths with total net flux equal to zero. Ephemeral re- 
gions emerge preferentially near supergranular boundaries, without influencing the 
flow pattern. After emergence, they most often fragment into se veral concentra- 



tions where their flux quickly merges into the pre-existing network. ISchriiver et al 



1997 ) give the average flux of an ephemeral region to be ~ 1.3x 10^^ Mx. iMartin 



( 1990f ) determined that the network concentrations are produced from the resid- 



uals of other flux concentrations. It is estimated that 90% or more of their flux 
originates from ephemeral regions, with the remaining 10% being contributed by 
intranetwork fields. Intranetwork (IN) fields are the smallest of the three types of 
small-scale magnetic flux concentrations. They typically last 2 to 4 hours. Their 
mean diameter is about 2 000 km, while their total flux content (~ 10^^ Mx) tends 
to be close to the detection limit of present instruments. They emerge within su- 
pergranule cells due to the emergence of small bipolar pairs. They are then swept 
towar ds the cell bo undaries b y radial flows, where they interact with network 
flelds (|Martinl[l98i) . Recently. [Centeno et al.l (|2007l) using observations obtained 
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with Solar Optical Telescope Spectropolarimeter (SOT/SP) on board the Hin- 
ode spacecraft presented clear evidence of an emerging IN magnetic loop showing 
strong horizontal magnetic signal and traces of vertical opposite polarities on each 
side of it. The lower parts of the loop were advected into the intergranular lanes, 
where they a ggregate to other magnetic field concentrations resulting in larger 
flux elements. Lites et aP (j2008l) based on the high spatial resolution and polari- 
metric sensitivity of the SOT/SP reported the existence of horizontal magnetic 
fields in the IN regions whose mean field strength surpasses considerably that of 
the vertical component. The vertical fields are concentrated in the intergranular 
lanes, whereas the stronger horizontal fields occur preferentially at the edges of the 
bright granules or inside them aside from the vertical fields. This kind of spatial 
relation between the horizontal and the vertical magnetic fields suggests that the 
IN field consists of small Q loops having their vertical roots in the intergranular 
lanes and connected by an horizontal field above granules. 

From the previous discussion it becomes clear that the magnetic flux patches 
within the magnetic carpet, while diffusing along the superganular boundaries 
due to supergranular flows, interact with one another by the following four main 
processes: Flux emergence which is the appearance in pairs or in clusters of new 
magnetic flux (with equal amounts of flux in both the positive and negative po- 
larities). Cancellation which occurs when concentrations of opposite polari ty come 
into contact and mutually lose flux ( Livi et aD. llQSSt iMartin et al.LTl985|) . If two 
or more concentrations of the same polarity merge together to form a single larger 
concentration, this is known as coalescence. Finally, fragmentation which is the split- 
ting of a large concentration of flux into several different smaller concentrations. 
Using a series of full-disc 96 minute MDI magnetograms, Hagenaar (2001) found 
that the time taken for all flux within the quiet Sun photosphere to be replaced 
("the flux replacement timescale") was around 14 hours. However, later studies 
determined this value to be an order of magnitude smaller, at just 1-2 hours 
(|Hagenaar et al.U2008[ ). 

The magnetic network formed by the driving of the magnetic flux to the bound- 
aries of supergranules is made up of clusters of kilogauss flux tubes. By definition 
the smallest observable structure of the magnetic fiux in the photosphere are mag- 
netic elements which have various field strengths and diameters close to or less than 
the best currently achievable spatial resolution (~ 150 km). Magnetic elements are 
very important because they not only structure the photosphere but also provide 
a link to the chromosphere and corona along which energy can be transported 
from the Sun 's interior to the outer laye rs. They were first observed as "mag- 



netic knots" (|Beckers and Schroteil Il968f) - small dark structures equivalent to 
what we now call "micr opores" (iTopka et al.L 1997|) - and as "filigree" , elongated 
chains of bright points ( Dunn and ZirkeiT 1973ft which have been resolved into 
strings of adjacent bright points by Mehltretteij ( 1974) . When observed at disc 
center they are usually best seen in the "G-band filtergrams" , which are obtained 
with a wide-band filter centered around 430.5 nm, named G-band by Fraunhofer. 
Th e first observations o f the n etwork at this wavelength window were carried out 
bv iMuller and Roudied (ll984D. The name "network bright points" (NBPs) was 
introduced fetenflo an d HarvevL [l985 iMulleil Il985|) for the bright grains which 
are visible in filtergrams and are arranged in roughly cellular patterns making up 
the network boundaries. Since bright points have been presumed to be the vis- 
ible manifestation of thin magnetic fiux tubes and thus related to the magnetic 
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Fig. 1 A SOHO/MDI magnctogram. The boundaries of supergranule cells are marked in 
yellow, and arrows indicate the supergranular flow pattern (Credits SoHO/MDI). 



elements, extensive work has been initiated and it is now well established that 
NBPs are brightness manifestations of the small, strong field magnetic elements 
that make up the magnetic network ( Muller and Roudieil , Il984t iBerger and Titlel 
Il996l I2OOII; iRouppe van der Voort et al.L 120051) . 

In stronger field areas, on the other hand, the increased number density of mag- 
netic flux tubes gives rise to large chromospheric areas called plages which con- 
sist of many bright elements (sometimes referred to as facular granules) packed 
together. Granules near network bright points and in plage regions are smaller, 
display slow er temporal evol ution and have a lower contrast than those found in 
quiet areas ( Title et al.l 19921). Th e granular pattern in these areas is referred too 
as "abnormal" ( Dunn et al. . 19741) . Flux concentrations take on structure, behave 
less passively and seem to affect the granular flow. The larger magnetic fllling 
factor in plage regions gives rise to large amorphous, elongated "ribbons" that 
contain high net flux. They are characterized by large flux concentrations having 
a brighter edge and an intensity depression in the center. The majority of them are 
not resolved into flux tubes, but instead show a substructure which is interpreted 
as an indication of a range of magnetic fleld strengths within the structure. Circu- 
lar manifestations of spread out ribbon structures have been called "flower s" . The 
"flowers" and "ribbons" are small flux concentrations ( Berger et al.L [20041) . while 
larger flux concentrations are called (micro)-pores and have a distinct dark center. 
There is a continuous transition between these features evolving from one to an- 
other depending on the amount of the magnetic flux that is collected or dispersed. 
The merging and splitting of "flux sheets" , which are concentrations of the mag- 
netic field due to the granular fiow, an d the transitions between the ribbons, fiow- 
ers, and micropores are described by iRouppe van der Voort et al. I (|2005|) . These 
authors found also that the plasma is basically at rest in these structures with 
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small concentrations of weak up-flow sites, while narrow sheets with down-drafts 
are found right at the edges of the magnetic field concentrations. 

Due to the continual motion of mixed polarity concentrations in the magnetic 
carpet and the different flux evolution processes, the quiet Sun photosphere is 
highly dynamic. Since magnetic flelds from the magnetic carpet extend up into 
the solar chromosphere and lower corona, it is expected that the quiet Sun chro- 
mosphere and corona are also highly dynamic responding to the movement and 
evolution of the different flux concentrations. Even more, the spatial distribution 
and rates of the different processes of flux evolution play a major role in deter- 
mining the distribution of heating events in the upper solar layers. The driving 
of foot-points of the small-scale magnetic fields on the surface of the Sun plays a 
fundamental role in sustaining the chromospheric network and leads to the build 
up of currents and/or the generation of waves, mechanisms considered as potential 
sources of energy for solving the coronal heating problem. On the other hand, the 
larger magnetic fill ing factor in plage regions g i ves ri se to more extended struc- 



tures. As stated bv lRouppe van der Voort et al. ( 20051) . the merging and splitting 



of fiux sheets and the continuous transition between ribbons and micropores, are 
difficult to reconcile with models that regard these structures as being composed 
of discrete fiux tubes that keep their identity over a long period of time. They 
also suggested that it should be more appropriate to characterize the temporal 
behaviour of the magnetic field in these areas as fiuid-like, where the "magnetic 
fiuid" is in continuous interaction with the field-free plasma, forming the extended 
concentrations. Due to the differences between quiet and active regions at the 
photospheric level described above, differences between the different underlying 
structures related to them should be expected. 

ft should be made clear at this point that the aim of this section is not to give 
an extensive review on the solar photosphere and its magnetic structure, but rather 
to provide the background for the overlying chromosphere and its structures which 
is the main topic of this review. The interested reader can be refer red to recent, 
more detai led revi e ws on the solar photosphere, like e.g. those by ide Wiin et al.l 
(|2009h and lMullej (|201ll) . 



3 The solar chromosphere 

The photospheric magnetic fields of different forms and strengths extend through 
the chromosphere into the corona providing the link between the different solar 
layers. As stated in the previous section in the photospheric layers the supergranu- 
lar flows carry the magnetic flelds to the edges of the supergranular cells resulting 
in large magnetic flux concentrations. The magnetic flux tubes that constitute the 
magnetic network expand upwards and appear as bright patches that constitute 
the chromospheric network. A strong spatial coincidence exists between the pho- 
tospheric magnetic network concentrations and the overlying network boundaries 
in the chromosphere. In Can H&K images, which sample the low chromosphere, 
NBPs show up with a larger bright ness enhancement over the surrounding area 
than they do in G-band (cf Fig. 2 of i Lites et al.l (199^), but with lower sharpness 
due to strong resonance scattering and possibly also due to increasing flux tube 
width with height (Fig. [2l upper row). Also high resolution images taken in the 
blue wing of Ha were found to display strikingly intense small-scale brightenings 
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Fig. 2 Images of a quiet-Sun area containing a quiescent filament obtained by the Dutch 
Open Telescope (DOT) in La Palma (Canary Islands). Top left: G-band. Top right: Can H. 
Bottom left: blue Ha wing at AX = —0.08 nm from line center. Bottom right: Ha line center. 
The left-hand one samples quiet net work, the right-hand o ne more enhanced network (figure 
courtesy of Dr J. Leenaarts, see also lLeenaarts et al.l ll2006l) ) 



at the locations where magnetic elements form the NBPs (Fig. [2] lower row, left). 
iLeenaarts et al.l ( 20061) showed the potential of these Ha bright points as proxies 
of the magnetic elements. 

The solar chromosphere being only about 2 500 km thick (~3") actually cov- 
ers the layer where the solar atmosphere turns from the gas-dominated pho- 
tosphere/lower chromosphere into the magnetic field dominated upper chromo- 
sphere/corona (for recent reviews on the solar chromosphere, a s well a s to it s 
coupling with the other solar layers the reader can be referred to iRuttenI (|2007l) . 
iRutten (2010), Wedemeyer-Bohm et al. (2009)). Chromospheric plasma, especially 
when seen in the center of the Hq line, is very clearly organized along fine-scale, 
dark, elongated structures that span out from regions of enhanced magnetic flux. 
Some of them are connected to neighbouring regions and others fading in between, 
while covering the IN cells (Fig. [2l lower row, right). Such filamentary structures 
portray the magnetic field topology and are consistent with the increasing dom- 
inance of the magnetic pressure over the gas pressure with increasing heights in 
the atmosphere. This dominance has as a result the upward expansion of the mag- 
netic fields which fill the space in the chromosphere and affect its dynamics. The 
relative importance of the pressures is given by the ratio of plasma pressure to 
magnetic pressure (also referred to as the plasma-/?) and is related to the forma- 
tion of the magnetic canopy zone. This zone, in which /3 = 1 and the sound speed 
equals the Alfven speed, separates the solar atmosphere into magnetic and non- 
magnetic regions and plays an important role in the refraction and reflection of 
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the upward propagating waves, as well as in converting ac oustic waves into other 
modes, such as fast a n d slow magnetoacoustic way esfse e e.g. Rosenth al et all 200 
iBogdan etlli |2003|) ) . iKontogiannis et al.l ([loiO^) and iKonto giannis et al I (|201o" 
have shown the important role the dark Hq fine-scale magnetic structures play in 
the formation of the magnetic canopy (for details see Section (H]). As one goes from 
the Ha line center to the wings, the picture gradually changes (see lower left panel 
of Fig. [2|). It is, usually, suggested that in dark structures observed in the Ha line 
wings the intensity is coming from lower layers than the core of the line, although 
it cannot be excluded that it might come from highly Doppler shifted features 
that are located higher in the atmosphere than the features observed in line cen- 
ter. The dark structures endings in the Hq blue wing image (Fig. [2] lower row, 
left) - which are part of the elongated mottles seen in Ha line center- are likely 
to emanate from photospheric bright points which are concentrated into clusters. 
The relationship between the ends of the dark structures and the photospheric 
bright points is presented in Section O 

From the above it becomes clear that the photosphere and chromosphere are 
coupled via magnetic fields and, furthermore, any description of the solar atmo- 
sphere in terms of average properties is incomplete without the consideration of 
the fine-scale structures and their dynamics. 



4 Fine-scale, jet-like chromospheric structures 

Mottles, spicules, fibrils and dynamic fibrils, straws, rapid blue excursions (RBEs) 
are some of the jet-like, fine-scale features that dominate the dynamic and highly- 
structured chromosphere. Direct, narrow-band or broad-band filtergrams have pro- 
vided the means for obtaining parameters such as lifetime, size, spatial distribu- 
tion, inclination, etc., of these structures. On the other hand, spectroscopic ob- 
servations in different lines are essential to extend our knowledge on the physical 
conditions, particularly the temperature, density and fiows. Determination of the 
physical conditions in the observed features is also based on the study of the spatial 
and temporal variation of the observed radiation intensity in different wavelengths, 
i.e. line intensity or contrast profiles. On the basis of observational findings various 
theoretical models and numerical simulations have been developed to explain the 
mechanism responsible for their formation. However, in spite of the remarkable 
advances made by the high resolution observations and by the development of 
theories and numerical simulations, their interrelationship, the determination of 
their physical parameters, the definition of their formation mechanism, as well as 
their possible role in the heating of the solar corona still remain uncertain. Am- 
biguities are mainly due to the differences in their appearance when observed in 
different spectral lines, but also at different wavelengths within the same line, to 
the different instruments used for their observation (which usually have different 
spatial and temporal resolution) , to the different methods used to infer their phys- 
ical parameters, especially velocities, to their relatively short duration and rapid 
changes and, most important, to the different names used by different authors to 
describe the same structure. 

Below we present a description of the appearance and the physical properties 
of these structures. We follow the traditional terminology, as well as the recent one 
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Fig. 3 A sample DOT Call H image obtained on November 4, 2003 showing numerous jet- 
like structures (spicules, active region fibrils, superpenumbral fibrils) clearly visible on the 
limb in addition to a large surge. The dark elongated structures near the limb are sunspots. 
At the bottom of the image thin bright structures, called straws, are emanating, from the 
chromospheric network (which is hardly visible in this image) , while around th e activ e regions 
several dynamic fibrils and penumbral fibrils are visible (from iTziotziou et al.l l|2005l) ) 



which has resulted from high resolution observations of some "new" , "spicular-like" 
structures. 



4.1 Off-hmb structures 



4-1.1 Spicules 



The traditional term "spicule" refers to the relatively thin, elongated jet-like struc- 
tures, seen at the limb of the Sun which appear as bright jet-like feature s against 
the d ark background of the solar corona. The term was first used by [Roberts! 
(|1945[) . ,0 describe jets which dominate the chromosphere seen at the limb shoot- 
ing up to 10 000 km heights. Spicule-like features, however, were first observed and 
descr ibed by Father Ang elo Secchi of the Vatican Observatory in 1877 (S ecchil 
1877| ). As first reported bv lLippincottI ( 19571) , spicules seem to show a group behav- 
ior, the so-called "porcupine" and "wheat" field patterns. In the former, spicules 
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seem to radiate outward from a common centre forming a 18 Mm wide configura- 
tion along the solar limb, whereas in the latter all spicules seem to have the same 
orientation over a width of 140 Mm. Initially, spicules were associated with quiet 
Sun regions at the limb, nowadays, however, it is believed that on-disc quiet Sun 
structures called mottles constitute only a fraction of the observed off-limb spicules 
and that active region fibrils, as well as superpenubral fibrils are seen as spicules 
when these regions cross the limb (see Fig. [3]). Despite their discovery over 130 
years ago, spicules are among the least understood phenomena of the solar chro- 
mosphere. Until recently, this was due to the lack of high-resolution observations 
(their diameters of a few hundred kilometers and lifetimes of a few minutes were 
close to the observational limits), as well as to the highly simplified and poorly 
constrained models used for their modelling. Several recent advances not only on 
the observational side, but also on the modelling side have provided significant 
progress on the understanding of these structures. However, we still lack definitive 
answers on several of their physical properties and interrelationships with other 
limb and o n-disk structures, wh ile models still face difficulties in describing them 
accurately. ISeckersI (|l968Lll972l) give comprehensive reviews of early work on mot- 
tles and sp icules, while a m ore recent review on the modelling of these structures 
is given bv lSterlind (|2000l) . 

The morphology and properties of spicules have been under study for the last 
sixty years although with rather limited interest after the early years of their 
characterization. Only recently, the interest of the solar physicists on them has 
been revived by the fact that high resolution instruments are being used. High 
resolution movies in the Ca ii H line obtained with the Broadband Filter Imager 
(BFI) of SOT on the seeing-free Hinode spacecraft have revealed a remarkably 
dynamic, spicule-dominated limb. Based on some physical pr operties derived from 
SOT/Hinode Can H observations. IPe Pontieu et al. ( 2007b[) suggested that there 
exist two types of spicules: a) "Type I" spicules, similar to the traditional spicules, 
which appear to rise up from the limb and then fall back, to be long-lived (3 - 
7min) and to exhibit longitudinal motions of the order of 20kms~^ and b) "Type 
H" spicules which seem to have much shorter lifetimes (50 - 100 s), higher velocities 

100kms~^), to be considerably taller and to exhibit only upward motion and 
then rapidly disappear (Fig. d]) . 

In that work the authors following the time evolution of the intensity along a 
fixed cut almost perpendicular to the limb (see Fig. they were able to show 
at least one clear example of a "Type I" spicule which shows a distinctive differ- 
ent profile than the steep linear streaks caused by "Type 11" spicules. However, 
since spicules are so dynamic and may move in and out of the fixed cut, such 
intensity profiles that cover the full time series cannot be used to reliably mea- 
sure the properties of spicules. Thus the authors had to rely on spicules that did 
not move trans versely so that their full life time could be covered by a fixed cut. 
Quite recently, Zhang et al.l ( 20121 ) re-analyzing the same data sets adopted by 
iDe Pontieu et all (|2007bi r identified and traced 105 and 102 spicules in quiet Sun 
and coronal hole (CH) , respectively. They claimed that they could not find a single 
convincing example of "Type H" spicules. Furthermore, they reported that more 
than 60% of the identified spicules in each region showed a complete cycle, i.e., 
the majority are "Type I" spicules. Thirdly, the lifetimes of the spicules in the 
quiet Sun and CH are 148 s and 112 s, respectively, but there is no fundamental 
lifetime difference between the spicules in the quiet Sun and CH reported earlier. 
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Fig. 4 The top panel shows spicules at the quiet Sun limb obtained with the SOT/Hinode 
in the Call H passband. The bottom panels show space-time (X - T) plots along the location 
indicated by a dashed line in the top panels, for both the original data and time difference 
data. The space-time plot is dominated by short-li ved vertical stripes (Type If spicules) and 
longer-lived parabolic paths (Type I spicules) (from lDe Pontieu et al.l ll2007bf )) 



In the same paper of lDe Pontieu et al.l (|2007b|) . the authors in order to reveal the 



two different populations applied two different filters to the time series. The first, 
low-frequency, filter was a Gaussian centered on 3mHz with a width of 0.5 mHz 
and brought out the "Type I" spicules, i.e. the relatively slowly evolving features 
that typically move up and down during their lifetimes of order 3 - 7min. The 
second filter was designed as a combination of a high-pass filter (> 15 mHz) with 
a low-pass filter placed at a frequency which is 2 mHz lower than the Nyquist 
frequency of the time series studied. This filtering isolates the longer lived Type 
I spicules from the more dynamic. Type II, spicules which appear and disappear 
on very short time scales. It has not been checked, however, if there exist spicules 
having lifetimes in-between the above mentioned timescales. 



iTavabi et al.l ( 201ll) have recently shown that there exist four types of spicules 



based on the spicule's diameter, ranging from 0.3" (220 km), 0.5" (360 km), 0.75" 
(550 km) to 1.15" (850 km). Unfortunately, detailed comparisons between the 
"classical" spicules and those called Type I and Type II is not easy to perform 
for several reasons among which the most important is the lack of spatial and 
temporal resolution of classical observations. Below we summarize properties of 
structures seen at the limb. 
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Spicule abundances In contrast to activity phenomena, spicules are omnipresent 
at the solar limb. Their abundance has been measured, as a variation of height 
and solar latitude. A point to emphasize is that the number of counted spicules 
depends very much on the quality of the observations. Thus their number, in total, 
is estima ted to reach ^ 10^ at zero height, while it is possible that this number 
i s larg er ( Beckersl 1 19721) . Most of the observational studies reviewed by Beckers 



(|1968n show that the spicule number increases with height reaching a maximum 



at a height and then decreases again. This is a result of spicules overlapping one 
another at lower heights and only few of them reaching large heights. An estimation 
of the variation of spicules' abundance with height gives aroun d 56 O OP sp icules 
on the full solar disk, around 5 000 km from the surface. Athavi ( 19761) and IZirinI 



( 19881) give 60000 and 70 OOP respectively, fo r the t otal number of spicules on the 



whole Sun at any time, ljudge and CarlssonI ( 2010l) assuming Type II spicules to 



be randomly distributed along the boundaries of circular supergranules into 200 
small bushes with a common "root" , and each bush to have eight spicules, they 
obtained 1600 spicules per supergranule. With each spicule having a diameter of 
0.1 Mm, they obtained an area fiUing factor of 0.015 and a total of 2x10^ Type II 
spicules on the Sun at a ny time. This value is some 20 times larger than the value 
given by iBeckerd (|l972l) . while Type I spicules are not included. Type I classical 



spicules being shorter are very difficult to count. They obtained this result by 
comparing Monte Carlo radiative transfer calculations with observations obtained 
by the SOT/Hinode in the Can H line. However. Sekse et al. (2012) questioned 
this estimate and argue tha t this numbe r is to o high and would imply a much 
more spicules than observed. iMoore et al.l ( 201lh surmised that at any given time 



there are ~ 50 Type II spicules present per supergranule, or ~1 Type II spicule 
per 2x10^^^ cm^ of surface area in quiet regions and coronal holes. All given values 
depend, of course, on the threshold intensity adopted for the spicule visibility 
and in which height from the solar limb spicules are counted. In some studies, 
a variation of abundance as a function of solar latitude has been found. More 
precisely, the number of spicu les was found to decrease from poles to equator 



(|Lippincottl Il957l: lAtha^l Il959l) . 



Spicule inclination The inclination of spicules is undoubtedly reflecting the local 
magnetic field topology. It should be mentioned that what can be measured is 
their apparent tilt with respect to the local vertical, since what we actually see 
is the projection of a spicule into the plane normal to the LOS and this leads to 
an underestimate of its true inclination. Measurements are also biased towards 
smaller inclination values, since the most inclined spicules are either not raised 
high enough to be visible or it is difficult to discer n the s horter ones due to super- 
position effects. Old studies, reviewed bv iBeckeri ( 19681) report a mean variation 



of ^ 2 0° of spicules' inclination from the local vertical. iHeristchi and MouradianI 



( I992I) measured the apparent tilts to the vertical of 843 spicules and from the dis- 



tribution of the inclinations they found an average tilt of 29° . It should be noted 
that their distribution was referred to an axis with a tilt of 14° in the yz plane 
as they found that spicul es tend to lean towards the equator. Recent calculations 
bv lPasachoff et al.l (|2009 l) did not find such a large tilt toward the equator. They 
measured an average absolute inclination to the vertical of 23°. Their statistics, 
however, contained only 40 spicules. A hint on the inclination of spicules has been 
achieved recently by the efforts to determine the magnetic field and its inclination 
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in spicules (see Section r4.1.ip . iLopez Arista and Casinil ( 2005h found a good cor- 



relation betwe en the orientation of the magnetic field lines and Ha spicules (see 
Fig. I12p . while IrVuiillo Bueno et all ([2005) found that the best fit to the observed 
Stokes profile is obtained for a magnetic field inclination oiO = 37°. |Pasachoff et al.l 
dlooa) pointed also out that it is possible that spicules' tilt varies with the solar 
cycle. Since spicules follow the pattern of the magnetic field, such a variation is 
possible. Considerable differences in the inclinations of spicules at different heli- 
ographic latitudes have also been reported. At the polar regions, spicules seem 
to have less variable inclinations and to be more perpendicular to the surface, 
whereas they attain more and more large r inclinations towards the equator as 
the latitude decreases. .Bonnet et al.l (1980) also found that in Ly-Q photographs 



of the solar limb, spicules appeared to be symmetrically distr i buted arotind the 
polar axis. There are several older studies ( van de Hulst . ITqSSI : iLippincottl |l957l) 



that associate spicules' inclination with that of the coronal polar rays, but such 
detailed correspondence should be re-examined with the current high resolution 
observations. 

Spicule lifetime and birthrate A spicule's lifetime is defined as the average dura- 
tion of spicule visibility. There is a large dispersion in the reported values of 
the lifetime of spicul es. Most studi e s agree that the li f etime of spicules is be- 



tween 2 and 12min (Roberts' '1945; 'Rush and Roberts', 'l954; 'LiDDincottl Il957l: 
lAlissandrakis and M aoris. 1 971 ; Cook et al. , 1984; Ge orgakilas et al. . 19991) with 
5 min considered as a typi cal mean value. The sam e range of values has been con- 



firmed for Ha spicules bv IPasachoff et aP ( 20091) who found lifetimes between 3 



and 12 min, with a mean value of 7.1±2.3min. The same authors find shorter life- 
times in TRACE 1600 A counterparts of Ha spicules pointing out a possible impact 
of the lower resolution of TRACE observations. The highest resolution available 
(from SOT/Hinode in the Can H passband) have attributed lifetimes around 3 - 
7 min for the classical spicules and lifetimes of th e order of 10 - 150 s to the ex- 
tremely short-lived ones termed Type H (|De Pontieu et al. l l2007b[) . lB^ckerd(|l968l) 
gives as a spicule birthrate the number 330 s~^. 

Spicule heights and widths The height of a spicule is conventionally measured from 
its foot at the photospheric limb up to where the spicule becomes invisible. It 
should be noted that what is finally measured is the projected height of the struc- 
ture on the plane of the sky. The height of a spicule is not a well defined quantity, 
since there is not a sharply defined upper boundary and its foot-point, as well as 
its location in front or behind the limb cannot be easily established, while its de- 
termination depends amongst others on the seeing conditions, line of observation, 
exposure time, etc. Especially, in the Ha line there is a chromospheric limb, which 
inhibits the tracing of the spicule to the photosphere. Apart from these difficulties, 
studies of t he distribution of th e apparent heights of spicules shows good general 
agreement. 'Becker^ (|l968l Il972l) reviewed older studies based on Ha observations, 



which gave estimates of average heights between 6 500 and 9 500 km. Sin ce then 
spicules have bee n reported to reach up t o 15 000 km (jCook et aL . 198^ in Ha 



but also in He H (iGeorgakilas et all Il999( ) . A recent analvsis bv Pasachoff e t al.l 
([2009) based on high resolution Ha observations of the solar North-Western limb 
combined with co-temporal TRACE images in the 1600 A channel gives appar- 
ent heights in the range 4 200 - 12 200 km with a mean of 7200±2 000km and 
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agree very well with the old results given bv iLippincottI ( 1957|) . They also identi- 



fied spicules common in both Ha and 1600 A image data sets and found that in 
1600 A images, they were taller (5 600 - 14 700 km) and correlate very well with the 
height in Ha. They concluded that the heights reached by spicules observed in the 
TRACE UV continuum are ~ 2 800 km greater than that reached in Ha. Similarly, 
spicular structure up to ^ 18 Mm has been reported obtained by TRA CE in Ly-a 
a n TRACE 1600 A at the n orth solar hmb (* Alissandra kis et al.l l2005^. According 
to lDe Pontieu et al.l (|2007b[) . the maximum lengths of spicules observed in Ca ii by 



Hinode vary from a few hundred kilometers to 10 000 kmwith most below 5 000 km. 
Type n spicules attain heights between 1 000 and 7 000 km and they are tallest in 
coronal holes. 

The diameter of spicules is probably the property most influenced by observing 
conditions, such as atmospheric seeing and the diffraction limit of the instrument, 
as well as overlapping effects that make it difficult to separate individual spicules. 
Older studies based on Ha and Can H and K lines estima te the spicules' diameters 
to range between 400 and 2 500 km ( Becker j . 19681 1972 ) , with most spicul e s hav - 



1960t iLvnch et al.i Il973l) 



ing diameters between 400 and 1 500 km (jPunn et al. 
which were characterized, at the time, by Beckers as the "most reliable" values. 
Recent high resolution Ha observations provide a more strict estimate , in the range 



300-1 100 km with a mean diameter of 660 km ( Pasachoff et al.l [2009!) which is in 



very good agreement with results by'Nishikawa (1988) who found an average diam- 
eter 615±250km. Pasachoff et al. (2009) also found that the width in the 1600 A 
TRACE channel is greater than the width measured with the Swedish 1-meter 
Solar Telescope (SST) at Ha by a factor of over 1.5, ranging between 700 and 
2 500 km. They argued that TRACE widths are impacted by the telescope's broad 
point-broad point-spread function, while the resolution of SST is approximately 
four times higher that that of TRACE. With the advent of high resolution in- 
struments, thinner structures are being observed (< 200 km) and some have been 
attributed to Type II spicules. They have also revealed t hat some Ca n spicules con- 
sist of finer threads (double thread structure) (S uematsu et al.ll2008l ). [Tavabi et al.l 



()201ll) give an extensive Table together with references (see their Table 1) of the 



various values obtained so far for the diameters of spicules. Furthermore, these 
authors based on Can H SOT/Hinode observations at the northern solar limb 
found that spicule diameters show a whole range of widths from 200 to more than 
1000 km, with distributions at four distinct diameters. 



Spicule temperatures and electron densities Spicule temperatures and densities can 
be estimated spectroscopically by using spectra taken simultaneously in several dif- 
ferent lines. It should be noted that temperature (and electron density) diagnostics 
of spicules from line intensities rely on data from constituents that are especially 
difficult to model. In order to interpret line profiles, non-LTE calculations have to 
be used which is not an easy task. Due to the difficulty of such calculations there 
are very few reports on temperature and density determinations in spicules and 
very large discrepancies on the reported values. Unfortunately, recent th eoretical 
results do not exist. Compilation of older theoretical results is given bv iBeckersI 
(|l972h . ] Beckers derived these parameters by assuming that a spicule is a cylinder 
situated vertically on the Sun and irradiated by the solar radiation having a di- 
ameter of 815 km. From the intensity of a spectral line he obtained a Te-Ne curve. 
He then compared the observed line intensities with the predictions of theoretical 
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Fig. 5 Relation between spicule temperature 
derived for different lines. Z stands for the He 
to hydrogen l ine intensity, R R for the ratio of 
spectra (from iBeckerj il972f) ) 



and electron density at different heights as 
10830 line, _R, for the ratio of 3888 A helium 
Call and He— D3 line as derived from eclipse 



models. Curves obtained by this way from different lines will intersect, giving a 
set of te mperature and density that is consistent with all observations (Fig. [5]). 
iBeckerd (|l972) derived the temperature from the observed intensities of Ha, Hen, 
and Can lines from simple non-LTE calculations. He reports a temperature of 
9 000K and an electron density of 1.6xl0^^cm~^ at a height of 2 000km and a 
temperature of ~ 16 000 K at 4 000 - 8 000 km, while for the electron densities he 
derived values 1.5xl0"cm~^ at 4 000km and 4.3xl0^°cm"^ at 8 000 km. 

iKrat and KratI ( 197l|) obtained simultaneous spectrograms of spicules in the 
lines Ha, H/3, D3 and Can H and K at heights varying from 5 000km to 9 000km 
above the limb, the slit being set parallel to the limb. They concluded that the 
emission in the different lines originates in different parts of spicules with different 
radial and turbulent velocities and that the electron density varies from lO^^cm^'^ 
to lO^^cm-^. 

lAlissandrakid ( 19731) analyzed simultaneous spectra of spicules at a height of 
5 400km above the limb in the Ha, H/9 and Can lines. Since the intensity of the 
Balmer lines is a very good indicator of the electron density, as they are almost 
insensitive to temperature, while the K line is sensitive to both the A'^e and Te, 
he combined the observations to obtain the temperature and density of spicules. 
He found that the values of the electron density for 36 spicules at a height of 
5400km range between 6xl0^"cm~'^ and 1.2xl0^^cm~^, with an average value of 
6xl0^''cm"^, while the range of the electron temperature was between 12 000K 
and 15 000 K with an average of 13 000 K. 

iKrall et al. I (|1976[) analyzed time sequences of simultaneous spectra of limb 
spicules in the Ha and Can H and K lines. Using the measured intensities in these 
lines they derived electron densities averaged over the entire visible lifetime of 
spicules of ~ 6xl0^°cm^'^ at the height of 5 000km and minimum and maximum 
values of ~ l.lxlO^cm"^ at 6000km and 2.xl0^°cm"^ at 10000km, respectively. 
They also found electron temperatures ranging between 12 000 K and 16 000 K, 
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while from profile half-widths, they concluded that turbulent velocities should be 
in the range between 12 and 22 km s~ ^ . 

iBraun and LindsevI (|l987l) used brightness limb observations at 100 and 200 /im 



and 2.6 mm. Use of far-infrared continuum observations are useful because the 
continuum emission is formed in LTE and the source function is just the Planck 
function, which in the Rayleigh-Jeans limit is proportional to the electron temper- 
ature. They calculated limb intensity profiles for a variety of spicule models and 
found that limb profiles are well fitted by spicules with an electron temperature of 
the order of 7 000 K up to heights of at least 7 000 km above the photosphere. Par- 
ticularly the 2.6 mm observations exclude spicule temperatures of 16 000 K below 
7000 km, because this temperature should create a substantial limb brightening, 
wh ich is not observed. 

iMatsuno and Hiravamal (jl988l) determined the height distribution of the ki- 
netic temperature of Ha spicules. The temperature was found to decrease from 
9 000 K at 2 200 km to 5 000 K at 3 250 km and to increase up to 8 200 K at 6 000 km. 
They suggested that the decreasing temperature with height might be related to 
the lateral expansion of rising spicular material along the expanding magnetic 
field lines, while the increasing temperature above 3 250 km may be due to heat- 
ing by the penetrating radiation in the Lyman continuum originating in the EUV 
em ission of the transition region-cor ona. 

ISocas-Nav arro and Elmora (|2005l) used multiline limb observations of spicules 
from different elements and, specifically, in the Can 8498 A and 8542 A lines using 
the Advanced Stokes Polarimeter (ASP) at the Dunn Solar Telescope (DST) and 
the He i multiplet at 10830 A recorded with the Spectro-Polarimeter from INfrared 
and Optical Regions (SPINOR), to derive some spicule properties. They found that 
the Doppler widths of the observed spectral lines are similar. Since the Ca atom 
is ten times heavier than the He atom, if the lines were broadened by microscopic 
thermal velocities, one would expect the He line to be more than three times 
broader than the Ca lines, which is not the case. Assuming non-thermal broadening 
they concluded that the electron temperature of the spicule should be lower than 
13 000 K to account for the observed line broadening. 

Evolution and velocities of limb structures For limb spicules, different and sometimes 
contradictory results on the velocities and their variations have been reported. One 
of the reasons is due to the different types of observations that have been used. 
Apparent velocities correspond to motions in the plane perpendicular to the LOS 
and could equally well result from the temporal variation of the excitation or 
ionization of the spicular material, or the propagation of a shock front. The gen- 
eral description based on measurements of apparent velocities is consistent with a 
spicule which rapidly elongates upwards with an average velocity of ~ 25kms~^ 
and reaches its maximum height within a minute or two after its initial appearance. 
The upward motion is usually fairly regular and continuous, the spicule rising and 
stopping abruptly at its maximum height. Subsequently, the spicule may either 
fade from visibility or else descend back to the low chromosphere with a velocity 
comparable to that of its initial ascent. Measurements of vertical pro per motions 
were made in the 1950's and 1960's by 'Rush and Roberts' ('1954') and 'LiDDincotlJ 
(1957), and have been re viewed bv [Becke rs (1968, 1972) and Bray and Loughhead 
( 1974f ). In an old work, iLiDDincottI ( 19571 ). after studying the apparent motions 



of spicules, reported that in half of them the ascending phase is followed by a 
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Fig. 6 Left. The apparent axial velocity distribution in spicules measured from the time 
variation of thei r heights. Right. The axial mass velocity distribution obtained from the Doppler 
velocities (from lHeristchi and MouradianI {1992)) 



descending one, while the other half is fading away; she was also the first to re- 
port that some spicules appear to rise from the same source several times. The 
major ity of descendin g spicules appear to retrace their original paths, but accord- 
ing to lBeckerj (Il972l) there are many which fall b ack along different paths, often 
forming an arc h. 'Heristchi and Mouradian' ('l99^ combining the observations of 
iRush and Robe rts (1954) and Lippincott (1957), which gave the variation of height 
with time of spicules, constructed the observed distribution of apparent velocities 
relating them to the angular distribution of spicules obtained from their axial ve- 
locity distribution (Fig.[6l left). They give as t he most probable velocity 20kms^^ 
and a range between Okms~^ and 60kms~^. IPe Pontieu et al.l ( 2007b|) from the 
analysis of chromospheric Ca ii H observations taken with the broadband filter on 
board the Hinode mission found that spicules sh oot upward at speeds between 20 
and 150 km s"^ Recentlv. [PasachofF et aP (|2009l) reported that in the majority of 
the 40 spicules examined, they found only ascending velocities because most of the 
spicules seemed to disappear at their maximum heights. They measured ascending 
velocities with a mean of 27.0±18.1 kms~^, a median of 25kms^^ and a range of 
3.0 to 75kms-i (Fig.B left). 

Apart from the velocities deduced from direct observations of the changes of 
the height of spicules with time, velocities can also be inferred from spectroscopic 
observations which give information on Doppler shifts. The Doppler shifts of the 
emission features can be interpreted as LOS mass motions of the spicular material. 
It should be noted that for a spicule inclined to the plane of observation, the 
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Fig. 7 Left. Histogram showing the distribution of apparent ascending velocities. Right. His- 
togram showing the distribution of the LOS velocities of th e tops of the spicules when they 
were at their maximum heights (from .Pasachoff ct alj ll200gf) ') 



component of the rise velocity corresponding to the sine of the angle of inclination 
would be observed. 



lAlissandrak"i3 (|l973l) using simultaneous spectra in Hq, H/3 and Can K at 



a height of 5 400 km above the limb obtained the LOS velocity distribution of 
50 spicules, as inferred from the position of the maximum of the line profile in 
the three lines. Thus he found an absolute average value of the LOS velocity of 
6.4k ms~^ for Hq, 5.6kms~^ for H/3 a nd 4.9 km s^^ for Can K. Using this tech- 
nique |Heristchi^ni^omadia3 ( I992I) measured the LOS velocity component of 



90 well-isolated spicules at 5 000-6 000 km above the limb. To improve the statis- 
tics they have also added some other available distributions (see their Table I). 
Relating the measured Doppler velocity distribution (which was approximately 
symmetric ranging between -30 to 30kms~^) with the angular distribution of in- 
clinations they obtained as a most probable axial velocity 40kms^^ (Fig.[6l right), 
i.e. ab out 20kms~^ la r ger th an the most probable apparent upward velocity. Re- 
cently [Pasachoffer^l] (|2009() measured LOS velocities from Doppler shifts of the 
tops of spicules found at their maximum heights (Fig. [3 right). The Doppler shifts 
were calculated by fitting Gaussians to 5— wavelength Ha spectra. They reported a 
mean absolute value of 5.1±5.1kms~^, a median of 3.8 km s^^ and a range of 0.3 to 
30.2 km s~^ This value agrees with the mean value of 6kms ^ quoted bv lBecker^ 
([1972). They also found that the distribution of the LOS velocities at the bases of 
the 40 spicules measured (where the base is defined as the location of the apparent 
limb) has a mean of 3.1±11.9kms~^ and a range of —17.8 to 51.2kms~^. As noted 
from their Doppler measurements, they regularly found oppositely directed mo- 
tions, i.e. bi-directional, which could be an indication of magnetic reconnection. 
One has to bear in mind, however, that whereas photographic observations give the 
vertical velocity of the spicules in the plane of the limb, spectroscopic observations 
give the velocity component along the LOS and, hence, in the direction normal to 
the plane of the limb. This fact can explain the discrepancy between the average 
apparent vertical velocity value of 25kms^^ and a n average value of ~_6kms~^ 
obtained for the LOS velocity of spicules. FoUowing lAthav and BessevI ( 1964[ ). one 



can show that the magnitude of the average Doppler velocity of spicules lying in 
a cone whose axis is normal to the LOS and whose apex angle is 26 (see Fig. [8]) is 
given by the formula: 
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Fig. 8 Resolution of the velocity V of a spicule along the LOS. Vjj is the LOS component, is 
the angle made by the spicule with the solar radius, and (p is the azimuth al angle of the pro jec- 
tion o f the spicule velocity vector onto a plane tangent to the limb ffrom lBrav and LoughheadI 

{Hzi)) 



vn = ^^^^^^^4^^^='-vsine (1) 

where is the azimuthal angle of the projection of the spicule velocity vector onto 
a plane tangent to the limb. At greater heights, the observed Doppler shifts will 
be positive or negative depending on whether the spicule is inclined away from 
or towards the observer. Thus for V=25kms^^ and ^=30° we get 8kms~^ for 
the average LOS velocity of spicules, which is close to the value reported from 

ob servations. 

iHeristchi and Mouradianl ( 1992fl . as stated above, found a large difference (-- 



20kms~^) between the maxima of the distributions of the axial velocities deduced 
from the measurements of the apparent and Doppler velocities of spicules. As they 
anticipated this difference is presumably due to the ionization of the hydrogen 
atoms, when the spicules, penetrating the corona, are heated. Thus if a partially 
ionized plasma parcel is moving along the axis of a spicule with a velocity of ~ 
40kms~^ and at the same time its front is being ionized and an ionization wave 
is moving backwards with a velocity of ~ 20kms^^ then the apparent velocity 
would be equal to 20kms^^. They also noted, however, that the differences in 
the obtained values may also be due to the fact that the data they used for their 
calculations were taken at different periods of the solar cycle, at different locations 
on the Sun, with different instruments and with different techniques. 

The height dependence of the LOS velocity in spicules has also been investi- 
gated by seve ral authors. The r esults obtained differ significantly from each other. 
For instance. iMouradianI ( 19651) has found that the Doppler shift decreases along 
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Fig. 9 SOT/Hinode Call H images showing spicules taken from left to right, on 2009, April 28 
at 16;24 UT, on April 29 at 02:24 UT and 03:28 UT. The over-plotted solid line indicates the 
position of th e SUMER slit, while th e dashed-dotted lines indicate the EIS raster field-of-view 
(FOV) (from lMadiarska et al.l J2011l'l ') 



the s picule, Beckers et al. (Il966l) found an increase of the Doppler shift with height, 



while fPasachoff et al. ( 1968h reported a weak or even no observable change of the 
shift with height. The main reason for these discrepancies are due to the fact that 
observations of spi cules at various height s simultaneously are related to instru- 
mental difficulties. iKulidzanish^ (|l980ll studied almost simultaneously height 
sequences of 69 spicules in the Ho line. He claimed that absolute values of LOS 
velocities increase linearly with height and that no variation of the sign of the LOS 
velocity along individual spicules is observed. 

It is quite evident that spicule motions and their time variations are essential 
in clarifying the mechanism responsible for their fo rmation. However , reported 
time variations of Doppler shifts are contradictory. iMouradianI (|l965h reported 
that the Doppler shift of spicules increases progressively with time, goes through 
a max imum and then decrea ses, but it does not change sign, being always up- 
wards. [Pasachoffer^l] ( 1968() found that Doppler velocities of a certain number of 
spicules do undergo quasi-periodic reversals. Recently, Wil helnJ (|200(]() , in a study 
of spicules observed in several EUV lines by SUMER on SOHO, re ported strong red 
and b lue shifts within a feature, which, furthermore, reverse sign. lMadiarska et al.l 
( 201l[ ) have recently analysed three large spicules (Fig. [9|) and found them to be 
comprised of numerous thin spicules which rise, rotate and descend simultane- 
ously forming a bush-like feature. Their rotation resembles the untwisting of a 
large flux rope. They show velocities ranging from 50 to 250 km s^^. These data 
will be further discussed in Part H of these review series. 

iNishikawal (|l988|) studied the proper motion of polar limb spicules using time 
sequence filtergrams of single wavelength bands centered on Ha-0.9 A or Hq-|-0.9 A 
(the passband was 0.5 A) obtained with the Domeless Solar Telescope at the Hida 
Observatory. After following successively from frame to frame four individual 
spicules he examined whether spicule motions can be interpreted as represent- 
ing ballistic or constant velocity. By fitting the height variation with a quadratic 
function of time (a constant acceleration model) he found large initial velocities of 
80 - 100 km s^^ and large acceleration of 0.45 - 0.65 km s~^. No disk features that 
can provide such initial velocities have been observed. He pointed out, however. 
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the possibility that such a source could be small enough that it is not observed 
with the current resolution and concluded that the ballistic model of spicule mo- 
tions cannot be rejected from the observational point of view. He also pointed out 
the possibility that a rising and a falling motion can be separately fitted with a 
constant velocity model. In this case, the constant velocities ranged from 30 to 
50 kms~^ 

IChristopoulou et al.l (|200ll) applied an image processing technique to high- 
resolution observations obtained with the DST of the Sacramento Peak Observa- 
tory. They estimated the proper motions of the apparent tops of several spicules, 
by performing a least-square fit to the height measurements as a function of time 
assuming a ballistic motion, i.e. that the plasma is subject only to gravitational 
forces. They obtained initial velocities larger than 40kms^^ in all cases examined, 
suggesting that the magnetic field or other forces should play an important role 
in the generation of spicules. It is clear from the above that the behaviour of the 
fiow along spicules remains an open question. 

Apart from the fiows along spicules, motions of spicules themselves parallel 
to the limb have been re ported. From sequence s of Hq spectra t aken with th e 
slit tangential to the l i mb.lPasachoff et al.l (|l968h . lNikolskvl (|l970j) . IWeartl (|l970l) . 
iNikolskv and Platoval ( 197ll) found evidence that some spicules change their po- 
sition along the slit relative to the average positions of other features thus ap- 
parently showing a real horizontal component of velocity. Nikolsky and Platova 
have concluded that the spicules execute quasi-periodic oscillations parallel to the 
limb with a peri od of ^ 1 min, characterist ic amplitude of l" and velocities about 
10 - 15kms^^. iKukhianidze et aP ( 2006h performed observations of solar limb 
spicules at 8 difTerent heights, i.e. ~ 3 800 - 8 700 km above the photosphere was 
covered. They found that ~ 20% of the measured height series showed a periodic 
spatial distributions in the Doppler velocities with a wavelength of 3.5 Mm and 
periods in the range of 35 - 70 s. They suggested that the spatial distribution was 
caused by transverse kink waves. Observations of such lateral or swaying motions 
have again been revealed recently from the unprecedented high spatial and tem- 
poral resolution images obtained with the SOT on-board the Hinode spacecraft. 
iDe Pontieu et al.l (|2007c') analyzed time series of Ca ii H images taken with a SOT 
broadband filter. They reported that many of the short-lived spicules undergo sub- 
stantial transverse displacements of the order of 500 to 1000 km during their short 
lifetimes of 10 to 300 s, while some longer-lived spicules undergo a transverse mo- 
tion with the displacement varying sinusoidally in time (Fig. llOp . They suggested 
that this behavior is strongly indicative of Alfven waves, where the term Alfven 
waves is used to describe incompressible transverse MHD waves that propagate 
along the magnetic field in an inhomogeneous medium. They pointed out, how- 
ever, that the observed waves could also be interpreted as MHD kink-mode waves, 
should a stable waveguide exist in the chromosphere. The fine structure and lateral 
motion of spicules has also been observed and reported bv lSuematsu et aP ( 2008h . 
These authors suggested that spicules can be driven by magnetic reconnection in 
unresolved spatial scale taking place at their foot-points. They also suggested that 
since most spicules emanate from a seemingly uni-polar magnetic region the rel- 
evant magnetic reconnection must take place in unresolved spatial scale contrary 
to the larger-scale jets associated with an emergence of a small bipole. 

Another unexplained feature that has been reported is that sometimes a "tilt" 
of the spectral line is observed in spectrograms, i.e., an inclination at a small angle 
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Fig. 10 (A) Spicules observed at the solar limb with SOT/Hinode in Can H. (B) A space- 
time plot (along the cut labeled 1 in (A)) of the Call H intensity. (D) A similar cut for the 
line labeled 2 in (A). The plot is dominated by a multitude of criss-crossed short linear tracks 
caused by spicular motion transverse to the magnetic field direction. Similar linear tracks are 
visible, as well as swings. The general characteristics (linear tracks and swings) of (B) and (D) 
are well reproduc ed by cuts that a r e gene rated from Monte Carlo simulations shown in (C) 
and (E) (from de lDe Pontieu et al.l ll2007cf) ') 



of the line profile of some spicul e s to the direc t ion o f the spectrograph's dispersion 
(|Beckersl , Il968l: iPasachoff etHI. Il968l: IWeartl. Il970t) . Such inclinations have been 
observed in the Ha and Ca H and K spectra of spicules and could be attributed to 
an actual difference in the LOS velocities between the two sides of the spicule. One 
obvious explanation would be that in these cases there is actually not one but two 
separate unresolved spicules with different Doppler velocities which are seen as one 
because of the low spatial resolution. Another explanation which was given is that 
this inclination may be interpreted as resulting from the rotation of the spicule 
about its axis which would produce differential mass motions. If, for example, the 
slit of the spectrograph is oriented across the width of a rotating spicule, with the 
side rotating towards the observer (towards the top of the slit) and the side rotating 
away (towards the bottom of the slit), then the spectral line will be inclined from 
the shorter w avelengths to longer wa velengths going from top to bottom along the 
spectral line. IPasachoff et al ] (|l968h gave an upper velocity limit for this spicule 
rotation at about 30kms~^. By considering a simplified model of a spicule as a 
column of gas 1000 km in diameter rotating with this peripheral speed they pointed 
out that it would have a velocity gradient of 60 m s^^ per km across it which would 
produce a spectral tilt of about 1°.9 in K and 2°. 6 in Ha. As they noticed this 
peripheral rotation velocity implies a centripetal acc eleration of ^ 1 .8km s~^, i.e. 
about 6 times greater than the value of solar gravitv. [AlissandrakisI ( 19731) . on the 
other hand, calculated a rotational velocity of ~ 8kms~^ and concluded that the 
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Fig. 11 Time series of sharpened images from Hinode Call H filtergrams. The cadence is 
5 s. The series shows that the s picule consisting of dou ble thread (indicated by white lines) is 
spinning as a whole body (from [siiematsu et alj ll2008l )^ 



rotation of spicules is almost negligible. Recently. IPe Pontieu et aD ( 20121) using 
high-quality observations from SST provided evidence that most, if not all, Type 
II spicules are characterized by the simultaneous action of three different types of 
motion: (1) field-aligned flows of order 50 - 100 kms^^, (2) swaying motions of 
order 15 - 20kms^^, and torsional motions of order 25 - 30kms^^. These last 
motions confirm the existence of twisting spicular motions. 

Hinode high resolution observations revealed th at most spicules show up a dou- 
ble thread structure during their evolution (Suema tsu et al-f l 20081) '). They also re- 
vealed that the separation of some of the double thread spicules changes with time, 
alternating between a single-thread phase and double-thread one (Fig. Ilip . This 
change in separation can be interpreted by a spinning of spicules as a whol e body 
(spin period of 1 - 1.5 min and velocity ~15kms~^). [Sterling et al.l (j2010f ) again 
from Can H SOT/Hinode observations reported that they found many spicules to 
expand laterally or split into two or more strands after being ejected. A possible 
explanation could be that many of the splitting or expanding spicules could be 
small-scale magnetic eruptions, analogous to coronal mass ejections. Consistent 
with this idea, the motion of the splitting spicules is similar to the spreading of 
the legs of filament eruptions. 



Magnetic fields in spicules The determination of the magnetic field vector in so- 
lar spicules is an important step to understand them. This determination can 
be achieved through specto-polarimetric measurements and theoretical modelling 
of the Zeeman and Hanle effects in suitably chosen spectral lines. Th e He i 10830 A 
triplet is a powerful diagnostic tool for chromospheric magnetic fields ( Truiillo Buenol . 
[2010 ). because its spectral signatures are sensitive to scattering polarization, while 
the Zeeman, Hanle and Paschen-Back effects make it a very useful magnetic indi- 
cator within a wide range of field strengths. 
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The Zeeman effect is the splitting of atomic energy levels (and of the associated 
spectral line) into several components in the presence of a magnetic field. In the 
linear regime this splitting is proportional to the magnetic field strength and the 
Lande factor of the energy level. The longitudinal Zeeman effect produces mea- 
surable Stokes- 1/ signals in the presence of weak magnetic fields, while to observe 
the signature of the transverse Zeeman effect on the Stokes Q and U profiles, field 
strengths of more than 100 G are needed. The Zeeman effect may produce polar- 
ized radiation that can be analyzed to infer the properties of the field. However, 
when the magnetic field is unresolved within the resolution element or cancels 
out when mixed polarities are present or too weak or absent then the splitting is 
negligible, there is no polarization signal. In these cases, the Zeeman effect is of 
no importance as a diagnostic tool of magnetic fields. There is, however, another 
physical mechanism, the so called Hanle effect, which allows us to "see" the "Sun's 
hidden magnetism" that the Zeeman effect is impossible to diagnose. The Hanle 
effect, contrary to the Zeeman effect, works in any topological complex weak-field 
scenario (even if the net magnetic fiux turns out to be exactly zero) . Atomic level 
polarization signals are produced by population imbalances and quantum coher- 
ence among the magnetic sublevels of the atom due to radiative transitions induced 
by the anisotropic incident radiation. The Hanle effect can be defined as any mod- 
ification of the linear polarization signals due to the presence of a magnetic field 
inclined with respect to the axis of symmetry of the radiation field. We should 
note, however, that for magnetic field strengths greater that 10 G (i.e. the sat- 
urated Hanle effect regime) , the Hanle effect is sensitive only to the orientation of 
the magnetic field vector, but not to its intensity. In this context, the Hanle and 
Zeeman effects can be suitably complemented for exploring magnetic fields in the 
solar atmosphere. 

The first direct empirica l demonstration of the exist ence of magnetized, spicu- 
lar material was achieved bv lTruiiho Bueno et al.l (|2005l) . They applied a combined 
Hanle-Zeeman diagnostic to the He I 10830A multiplet to spectropolarimetric ob- 
servations obtained with the Tenerife Infrared Polarimeter (TIP) at the German 
Vacuum Tower Telescope (VTT) on Tenerife, Spain. The spectrograph slit was 
placed off-limb parallel to the solar limb at an atmospheric height of « 2 000 km. 
They detected non-zero Stokes U profiles which according to the Hanle effect the- 
ory is the observational signature of the presence of a magnetic field inclined with 
respect to the local vertical. To correctly model the Stokes I profiles, the authors 
solved the radiative transfer equation assuming an optically thick atmosphere. 
They found from the analysis of quiet-Sun limb spicules at a height of 2 000 km 
above the photosphere that the best fit to the observed Stokes profiles is obtained 
for a magnetic field strength of B 10 G and an inclination angle 6 ^ 37°. The 
authors noted that the observed Stokes profiles are due to an averaging along the 
LOS and thus the possibility of stronger fields occupying a small fraction of the 
integration volume along the LOS cannot be excluded. 

ILopcz Ariste and Casinj ( 2005^ used spectropolarimetric observations of spicules 
in the He I D3 line done with the ASP at the DST of the Sacramento Peak Obser- 
vatory. The spectrograph slit was placed parallel to the solar limb at ~ 3 500 km 
above the Ha limb. The data consisted of full Stokes profiles showing significant 
broadening due to non— thermal effects. They constructed a database of synthetic 
profiles with different magnetic field parameters. They then applied the Principal 
Component Analysis (PCA) inversion to a series of spectropolarimetric profiles 
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Fig. 12 In the background, is the slitjaw image in Ha showing the position of the slit (horizon- 
tal) over the spicules. Superimposed on the same column (with a vertical offset) the inverted 
solution for the magnetic field vector for each point on the slit is plotted. The direction of the 
lines is given by the azimuth in the LOS reference frame, and their length by the field strength. 
The level of gray of each line is scaled according to the intensity of the observed feature in the 
slit-jaw image (from Lopez Ariste & Casini (2005)) 



which is performed by searching for the best match of the observed profiles within 
the precomputed database of profiles. They found a good correlation between the 
magnetic field orientation and the spicular structures visible in Ha (see Fig. I12p . 
However, due to a 90° ambiguity which they attributed to the small significance 
of Stokes V in their data, a solution giving a magnetic field approximately perpen- 
dicular to the axis of a spicule cannot be excluded a priori. They concluded that 
weak fields (of the order of 10 G) are more common in spicules, although stronger 
fields above 30 G should not be excluded in a minority of cases. We must note, 
however, that these sp icules emanated from an ac t ive pl age. 

In the same year, ISocas-Navarro and Elmord (|2005 l) obtained off-limb spec- 
tropolarimetric observations in the Hei 10830 A triplet using SPINOR in the Can 
8498A and 8542A lines using the ASP at the DST. They found that the ratio of the 
Q and U signals is different for the Ca and He lines implying a different orientation 
of the polarization plane. They concluded that this result provides clear evidence 
of the presence of magnetic fields in spicules by the following reasoning: if there 
was no magnetic field, one would observe some degree of scattering polarization. 
The plane of polarization of this signal is given by the geometry of the scatter- 
ing process and would be the same for all three lines observed. If, however, the 
scattering atoms are embedded in a magnetic field, then the Hanle effect would 
rotate the plane of polarization by an amount that depends on some parameters 
of the transition. In this case, the orientation of the linear polarization would be 
different in d ifferent lines, a s was observed. 

Recently. I Centeno et al.l (|2010l) carried out spectropolarimetric observations of 
off-limb spicules at various distances from the visible limb in the He I 10830A 
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triplet with the TIP at the VTT on Tenerife. With these particular observations, 
they detected clear Zeeman-induced Stokes V profiles. Such profiles are due to a 
net LOS component of the magnetic field. T hey used the inv ersi on cod e HAZEL 
(HAnle and ZEeman Light) developed bv lAsensio Ramos et all ( 20081) to invert 
the Stokes profiles caused by the joint action of atomic level polarization and the 
Zeeman and Hanle effects in order to infer a number of atmospheric properties 
together with the magnetic field vector. They inverted the Stokes profiles for all 
the positions along the slit for two data-sets taken one hour apart from each other. 
In a small region of the slit they found magnetic field strengths as large as 48 G, 
where the Zeeman-induced Stokes V signals were particularly large. They also 
pointed out that the determination of the longitudinal component of the magnetic 
field depends on the distribution of the orientation of the magnetic fields which 
may produce a cancellation of the Zeeman-induced Stokes V signal when magnetic 
field vectors are not aligned, or a significant non-negligible Stokes V signature (and 
thus a significant Bj^qs component) when the magnetic field lines team up. 

Another method that has been used for the det ermination of magnetic fields i n 
spicules is based on the MHD coronal seismology ( Nakariakov and Ofmanl l200lh . 
Properties of waves and oscillations (e.g. amplitudes, periods), as well as physical 
parameters of the medium (e.g. temperature, density) determined from observa- 
tions can be connected through MHD wave theory to the unknown magnetic field 
strength and transport coefficients. 

Magnetic fields in spicules can also be inferred from observations of waves. Sev- 
eral observations have shown oscillatory transverse displacement of spicule axes. 
T his d isplace ment can be caused by two types of waves: kink or Alfven waves. 
In IZaq arashvil i and Erdelvil (|2009f) . there is a discussion on difficulties associated 
with the Alfven wave scenario. Most authors have interpreted this kind of oscilla- 
tion in spicules as the observational signature of propagating kink waves. It is well 
known that transverse kink waves in ffux tub es anchor e djn the photosphere c an be 
gener ated by buffeting of granular motions (jRobertsl il97si : [Hasan and Kalkofenl 
Il999l) . The propagation of kink waves can be traced either by direct observation 
of the tube displacement along the limb or spectroscopically by the Doppler shift 
of spectral lines when the velocity of the kink wave is polarized in the plane of 
observation. If spicules are considered as ffux tubes anchored in the photosphere 
then the observed transver se displacements of t heir a xes can be interpreted by the 
propagation of kink waves. ^Kukhiani dze et all ( 2006[ ) analyzed Hq spectra of limb 
spicules obtained at different heights. They found that ~ 20% of the measured 
heights showed a periodic spatial distribution in the Doppler velocities which they 
attributed to kink wave propagation. Wave periods were estimated as 35 - 70 s 
based on the expected kink speed in the chromosphere (50-100 kms^^). Estimated 
wavelengths at the photospheric level are comparable t o the spatial dimens i ons o f 
granular cells, suggesting a granular origin of the waves. ISingh and Dwivedil (|2007l) 
used these observations and applied the method of MHD seismology including the 
effects of gravitational stratifica tion. The y estimated magnetic fields in spicules 
in the range 8 - 16 G. Zaqar ashvili et al.l (|2007i) analyzed the same observational 
data, estimating the kink speed in the ran ge 90 — 115kms~^ they obtained mag- 
netic field strengths in the range 12 - 15 G iKim et al.l (j2008h used high resolution 
Can H limb observations obtained by the SOT/Hinode instrument to determine 
oscillation parameters of spicules, such as periods, amplitudes, transverse veloci- 
ties, wavelengths and wave speeds. They interpreted the observed oscillations as 
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MHD kink waves and adopting spicule densities in the range 2.2 10^^^ kg m^'^ 
- 4.0 io~io kg n i ~^ th ey estimated magnetic field strengths in the range 10 - 
76 G. IVerth et all ( 201lh used also SOT/Hinode Can H observations and magne- 



toseismology to determine the vertical gradient in both magnetic field and plasma 
density in a spicule by studying the change in velocity and phase speed with height. 
They found that the magnetic field decreases in strength by a factor of 245 be- 
tween the photosphere and low corona, while the plasma density decreases by a 
factor of ~ 1000. However, in all of these works, various assumptions have been 
used. 

The measured values for magnetic field and inclination, therefore, are lower 
than the values needed for efficient p-mode excitation in inclined flux tubes. A 
closer investigation of the mechanism, the assumptions made in the analysis cou- 
pled with a larger number of spicule magnetic field measurements are required in 
order to resolve this issue. 

4- 1-2 Oscillations and waves in spicules 

Oscillations and waves in spicules have been detected by both imaging and spec- 
troscopic observations. Reported observations on oscillatory and wave phenomena 
in spicules, as well as vi ews and discussions about their i nterpretation can be found 
in the recent review bv lZaaarashvili and Erdelvil (|2009l) . 



4.2 On-disk chromospheric structures 

The chromosphere is the layer where the plasma-/3, changes from above to below 
unity, signaling a shift from hydrodynamic to magnetic forces as the dominant 
agent in the structuring of the atmosphere. Due to the combined effects of the 
pivotal role played by the magnetic field and the small-scale gas thermodynamics, 
this part of the atmosphere is characterized by an impressive amount of fine-scale 
structures. In active regions, long and relatively stable fibrils are observed, together 
with shorter dynamic fibrils, while in quiet Sun regions dark mottles emanate from 
the magnetic network. Whether or not these structures are similar and are driven 
by the same mechanism, as well as their relationship to the limb spicules, has been 
the subject of long-standing discussions for several years. 

The line which has been widely used for chromospheric observations, especially 
following the development of the Lyot filter, is the Balmer Ha line. Thus despite 
the complicated formation mechanism of this line and the difficulties in the inter- 
pretation of Hq observations, much of what we know about the fine-scale on-disk 
structures (e.g. mottles, fibrils) has been obtained from observations in this line. 
Substantial work has also been done by using the unique capabilities of the Multi- 
channel Subtractive Do uble Pass (M SDP) spectrograph which can acquire spectra 
over an extended FOV (iMeinl. |200 . The Can H and K resonance lines, as well as 
the Can infrared triplet (IR) (8498 A, 8542 A and 8662 A), have also been used for 
the solar chromosphere diagnostics. The former are the broadest lines in the visible 
spectrum and sample a wide height of the solar atmosphere. On-disk filtergrams in 
the Ca II H and K line are usually obtained with bro ad-band filters and differ sig- 
nificantly in appearance from those obtained in Ha. iRuttenI (|2007|) has remarked 



that dark structures in Ca ii H and K filtergrams are barely visible. Observations 
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in these lines close to the limb show long thin emission features called "straws" . 
iReardon etaP (|2009|) . however, have shown that the lack of well defined chromo- 



spheric features in the Can K fU tergrams, such as fibr ils and mottles, is essentially 
due to observational limitations. iPietarila et aT 1 (|2009') found that bright fibrils are 
ubiquitous in the Ca ii K line provided the spatial resolution is sufficiently high and 
the wavelength band is sufficiently narrow to avoid contamination by photospheric 
radiation. The Can IR and Ha lines differ substantially in their sensitivity to the 
temperature. The former are, like Ha, subordinate, but whereas the lower level of 
Hq is coupled to the hydrogen ground level via the extremely strong Ly-a radiative 
transition, the lower level of the Can IR lines is metastable and only coupled to 
the Can ground level via electronic collisions. For more extensive disc ussion on the 
forma ti on characteristics of thes e important chromosph eric lines see ICauzzi et al.l 
( 2008h , ICauzzi et al.l ( 2009() and iLeenaarts et al.l ( 2012h. High resolut i on im aging 



spectroscopy at fast cadence in Can IR was initiated bv lVecchio et al ] (|2007l) with 
the Interferometric BIdimensional Spectrometer ( IBIS) installed at the DST of the 
National Solar Observatory (NSO). ICauzzi et al.l (j2008l) using the same instrument 
obtained monochromatic images at several wavelengths within the Can IR line, 
analyzing several structures. They found that the appearance of the structures 
in this line is strongly reminiscent of the structures seen in Ha. As in Ha, the 
structures observed in Can IR lines differ noticeably between active regions and 
the quiet Sun. 

4.2.1 Mottles 

The traditional term mottles refers to rapidly changing hair-like jets observed in 
quiet Sun regions on the solar disc, usually in the Ha or the Can lines. They have 
long been recognized to be one of the basic elements that constitute the inhomo- 
geneous chromosphere. Several authors tend to refer to these features as spicules, 
although spicules do not always necessarily represent only quiet Sun structures 
seen at the limb. Mottles are organized in a complex geometric pattern over the 
solar disk outlining the boundaries of the chromospheric network which is more 
prominent in the Ca lines. They cluster into: a) either small groups called chains, 
consisting of almost parallel structures emanating along the common boundary of 
two supergranular cells or b) in larger groups called rosettes, consisting of usually 
radially expanding structures around the common boundary area of three or more 
supergranular cells. Tanaka (1974) found about 30% of all the dark mottles are 
double or can be resolved into double structures when viewed with high resolution 
Ha wing observations. 

Mottles appear as dark (absorbing) against the disk when observed in the wings 
of Ha and especially enhanced in contrast at Ha±0.5 A. Near the centre of the Ha 
line they are much less distinct. Very often in the past it has been reported that 
at the same locations dark as well as bright mottles are observed. However, even 
such basic considerations as whether dark and bright mottles are the same feature 
seen at different heights or not, remains unanswered. So me authors suggested tha t 
a bright mottle is the base of a dark elongated one (jBanos and Macrisl. Il970l) . 



whereas others claimed that bright and d ark mottles are distinctly different phe- 
nomena ( Alissandrakis and Ma cri3 . Il97lh . Today with the achievement of higher 
resolution observations one has the impression that what was called bright mottles 
is simply the bright background below the dark mottles. It is very important to 
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Fi g. 13 Hq i n tensitie s of the mottles vs the logarithm of the electron density A^e as computed 
by IGiovanellil lll967bl ) (in cgs units). Different curves correspond to different temperatures. 
Geometrical thickness of a mottle is D = 2 Mm 



note that the resolution of observations affect drastically the observed morphology 
of the chromospheric fine structures. Since mottles are usually grouped in rosettes 
or chains one will fail to discriminate between different structures when using low 
resolution observations, since each thread and its separation are a few tenths of 
arcsec wide. With the increasing capability of instruments such as SOUP, CRISP, 
IBIS or ROSA, we are able to detect fine structure with dimensions at the limit 
of the resolution. 



Lengths and widths of mottles iBrav and LoughheadI (|l974 in their monograph, 
summarize their properties. Their horizontal di mension ranges between 725 - 
10 000 km and they exhibit a variety of shapes. ISawved ( 19721) categorized the 
shapes of mottles into oval, round, filament, arch and lumpy and gave horizontal 
dimensions between 3"-10.6" (2 000 - 7500 km). Although generally elongated, 
the larger mottles tend to have more irregular shapes. Observations taken towards 
the limb, at large heliographic latitudes show mottles reduced in size, with widths 
of the order of l" and dominantly oriented towards the limb. 



Lifetime of mottles The lifetime of mottles is difficult to determine. Although the 
general configuration of mottles around bright points may be preserved for more 
than 15min, features change within minutes. Old studies give values ranging be- 
tween 3 and 15min, but differences in re s olutio n may affect the measurements. 
A more recent analysis bv lBratsolis et aL 1 (|l993l) determines the mean lifetime of 
mottles between 13-14 min. 
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Spectroscopic properties of mottles iGiovanellil (|l967allbh did pioneering extensive 
NLTE computations of the Ha line contrast profiles relevant to chromospheric 
fine-scale modelling. In his computations he assumed ID slabs illuminated by the 
surrounding atmosphere, having a geometrical thickness 2 000 km and uniform Te 
andiVe. For mottles' diagnostics, an important figure is Fig. II 31 given bv lGiovanellil 
(|l967bl) . In this figure, individual curves represent the Ha line-center intensities as 
a function of electron density for different temperatures. From this figure one can 
conclude that portions of the curves below the horizontal line correspond to dark 
mottles, while those above the line are structures considered as bright mottles. 
An important ambiguity one can notice is that for each negative contrast, for a 
given temperature, two values of electron d ensity are obtained. Using the observed 
contrast measurements, IGiovanellil (|l967b|) has concluded that: a) for dark mottles 
T < 10 000 K, A^e ~ 2xl0^^cm~^, b) for less opaque mottles T > 20 000 K, iVe ~ 
10"cm"3 and c) for bright mottles T < 25 000K, Ne > 5xl0^°-10"cm~3. 

iHeinzel and Schmiedeil ( 19941 ) used Ha line profile observations of bright and 
dark mottles together with a grid of prominence-like NLTE models of lGouttebroze et al.l 
(fl993l to derive the physical conditions in these structures. In these models, the 
structures are considered as vertically-standing ID slabs irradiated from both sides 
by an isotropic incident radiation. They have shown that higher-pressure models 
{pg « 0.5dyncm~^) with temperatures around 10"* K can explain the profiles of 
both dark and bright structures. However, it should be pointed out that the NLTE 
models they used are more suited for prominences located at a height of 10 000 km. 

lAlissandrakis et al 1 (|l990^ . [Tihopoula et al.l (|l993|) . and fTziotziou et d1(|2003h 
used Ha time series observations of mottles obtained by the MSDP and the cloud 
model (see Section [8.ip . They derived four adjustable parameters of the model, 
i.e. the source function, the Doppler width ziAo, the optical thickness tq, and the 
LOS cloud velocity v. Although the Ha observations were obtained in different 
times and regions the results they quoted are very similar. According to the last 
paper these structures are mostly optically thin with an optical thickness peak- 
ing at ~ 0.9. This is refiected on the source function distribution which peaks at 
~ 0.156 (the line center quiet-Sun Ha intensity relative to the continuum being 
0. 169) , since optically thin structures allow more of the background radiation to be 
transmitted through them. Their Doppler width distribution shows a peak around 
0.43 A and has a mean value of 0.35 A. The corresponding histograms of the four 
cloud parameters (including the LOS velocity which is described in the following 
sub-section) are given in Fig. 1141 

iTsiropoula and Schmieded (|l997l) proposed a method which allows the estima- 
tion of several physical parameters of dark structures once the four parameters 
described above have been derived through the cloud model for Ha observations. 
These parameters are: the population densities at levels 1, 2, 3 (A'^i, N2, N^) of 
the hydrogen atom, the total particle density of hydrogen Nj^, electron density 
A^e, electron temperature Te, gas pressure p, total column mass M, mass density 
p, and degree of hydrogen ionization xh- Mean values of these parameters, as well 
as their standard deviations, estimated from a set of several mottles are given in 
Table 2. The temperature has been calculated from the Doppler width obtained by 
the cloud model. They have assumed two different values for the microturbulent 
velocity, (i.e. lOkms"^ and 15kms~^), showing its effect on the values of the 
temperature and pressure. An assumed value of equal to lOkms"^ lead to 
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Fig. 14 Histograms of velocity v, optical thickness rp , Doppler wid t h zlA p and source function 
S derived from the cloud model for mottles ffrom pTziotziou et al.l ||2003| )1 

Table 1 Physical parameters of dark mottles from lTsiropoula and Schmiedej lll997l ) 



Parameter 




Average 
value 


Standard 
deviation 


A''i, cm~3 




1.6 10^" 


8.3 10'' 


N2, cm-3 




1.4 10"* 


1.1 10* 


N3, 




1.6 10^ 


1.3 10^ 


Nh cm~3 




5.1 10^° 


2.1 10^° 


Ne cm-3 




3.4 IQio 


1.5 10i° 


M, g cm~^ 




2.2 10~5 


9.4 lO"*^ 


p, g cm-3 




1.1 IQ-" 


4.7 10"" 


T, K (for it 


= 10 kms-i) 


1.4 10'' 


9.2 10^ 


T, K (for it 


= 15 kms^i) 


1.0 10" 


7.7 10^ 


p, dyn cm~'^ 


(for it = 10 kms-i) 


0.20 


0.1 


p, dyn cm~^ 


(for It = 15 kms-i) 


0.14 


0.1 


Xh 




0.65 


0.1 



1.7 times greater values of temperature and pressure than the values obtained by 
assuming equal to 15kms^^. 

Spatial and temporal evolution and velocities of dark mottles Among other physical 
conditions, the study of the flows along these structures is essential because it can 
add to the understanding of the mechanism driving them. However, the problem 
is difficult and the present situation is rather confusing because to determine a 
mottle's motion and its variation with time is a very delicate and difHcult matter. 
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This is not surprising in view of the difficulties of recovering information on the 
velocities from line profile observations formed from moving material embedded 
in an absorbing medium at rest. Stable and good seeing conditions are required, 
while wide wavelength coverage is necessary in order to follow the changes of the 
LOS velocity. Furthermore, careful co-alignment of images from different times, 
but also of different wavelengths is required. 

Another problem is that the apparent shape of mottles is not simple and is 
changing with time. Some are very straight, while others are curved. Some mot- 
tles are very thin and others are very thick. Some are tapered and oth e rs are 
not. In hig h resolution Ha line wing filtergrams, as described bv iTanakal ( 19741) 



and iDara et al., (.1998.), thev show a twined or multiple appearance. It is worth 



noting that Suematsu et al.l ( 2008f) found the similar twined structure for limb 



spicules (see Fig. Ilip . Furthermore, the mottles in the Ha blue- wing are thinner 
and closer to the netw ork boundaries than those observed in the Ha red-wing 
(jSuematsu et al.ll 19951) . 



A common practice to infer LOS velocities is to measure Doppler velocities 
by determining the wavelength of a point midway between two positions of equal 
intensity in the absorption line profile. In this technique the moving structure is 
assumed to produce its own profile which is simply shifted with respect to the 
background profile by an amount AXj = Xqv/c, where v is the LOS velocity. 
The determination of the velocities by this technique usually yields inconsistent 
results due to the lack of knowledge of the effective height of formation of the lines 
considered. Another technique used for the determination of LOS velocities, when 
intensities on either side of a line profile are known, is based on the red-blue wing 
subtraction 

I{+AX)-Ii-AX) 
I{+AX) + I{-AX) 

This technique gives the so-called "Doppler signal" . According to this formula, 
positive Doppler signal denotes upward motion of absorbing material. Of course, 
one must be aware that the derived values give only a parametric description of 
the actual velocity field, they can give, however, a qualitative picture of upward 
and downward moving material (Tsiropoula. 200Q) . 

Although mottles show different appearance and are someti mes displaced lat' 



erally among filtergrams of different wavelengths in the Ha line (jSuematsu et al.K 
[l995i) . we should note that no one has analyzed their horizontal motions. iBeckerg 
(|l963l) reported that dark mottles, after disappearing in the Ha-0.5 A image, often 
became visible in the Ha-|-0.5 A, consistent with up-fiows followed by down-fiows. 
iBhavilail (|l965l) has shown that some mottles are visible in the blue wing, while 
others appear only in the red wing. Upward and downward ve locities, as inferred 
from Hq spectra, have been reported bv , Suematsu et al.l ( 19951) . These authors af- 



ter careful co-alignment and examination of a time series of high resolution images 
at Ha-0.65A, Ha line center, and Ha+0.65A and Doppler signals proposed the 
following description for the evolution of disk mottles: at first they appear as faint 
dark features at the Ha blue-wing filtergrams, elongate and increase in contrast, 
and then fade out as a whole in this wing. A few minutes later, they appear again 
as dark features at the Ha line center, and show increased contrast and length. 
Finally, they appear at the Ha red-wing as faint, but long dark features, and 
shrink towards the footpoint, increasing their contrast and then fading out. They 
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Fig. 15 Time slice images of cloud velocity v (left) and temperature (right) along a dark 
mottle. On the velocity images the black contours denote upward velocities, white contours 
downward velocities, while the thick gray line represents the zero velocity contour (from 
iTziotziou et al.l ||2003I1 ^ 



concluded that Doppler signals confirm this pattern, i.e. upward radial velocity 
during the extending phase of the structure and downward during the receding 
phase or of an up— and— down moving jet as a whole. They also found that in some 
cases, the red-wing mottles appear near the foot-points of the blue-wing mottles 
in their extending phase. Cauz zi et ah. (,2009,) , on the other hand, compared profile 
minimum intensities and Doppler shifts of the minimum intensities of Ha profile 
observed with IBIS. Although there are large scale spatial similarities between the 
Hq profile-minimum map and the corresponding Doppler shift map, they found 
two dissimilarities: a) there is a lack of pixel-by-pixel similarity between the two 
maps and b) mottles in Doppler-shift maps tend to be shorter and are dark or 
bright (moving upward or downward) without obvious relation to the contrast in 
the profile-minimum image, except for the brightest and darkest features whose 
co-spatiality indicates exceptional motions in their lower ends. In the correspond- 
ing Can infrared line observations of intensities and Doppler shifts, they found 
slightly higher coherence but with negative correlation, i.e. higher intensity values 
at locations with la r ger bl ue-shifts (dark in the maps). 

ISuematsu et al.l ( 19951) examined also the proper motion of the apparent top 
of mottles. They showed that the tops of some mottles approximately follow a 
parabolic trajectory, when they are traced through the three wavelengths along 
the Ha line profile. Fitting the trajectory with a quadratic function, they derived 
an initial velocity and a deceleration of proper motion. Furthermore, they gave 
a method to derive a tilt angle from the vertical, assuming that mottles move 
on inclined straight flux tubes and are decelerated only by solar gravity. They 
found initial velocities ranging from 10 to 50 km s^^ with a mean of 28 km s^^ and 
decelerations ranging from -0.2 to -0.02 km s"^ with a mean of — 0.07 km s , 
which are much smaller than the solar gravitational acceleration of 0.274 km s^^ 
(|SuematsuLll998[) . 

If mottles follow a ballistic trajectory (material motion under a constant grav- 
itational field) , we can expect to find some relationship between their lengths and 
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lifetimes. For example, if the mottle moves along a rigid magnetic flux tube of 
a given inclination angle, with an initial velocity along the flux tube of vq, de- 
celerated by solar gravity g, we have a relation between the maximum length L 
projected on a horizontal plane and a half-lifetime T of 



Vo 



2 VT^^ (3) 

where Tp = vo/g is the lifetime of a vertical mottle whose apparent length is zero. 
This equation implies that longer lived mottles have larger ap parent lengths, if ev- 
ery mottle has a similar initial velocity. ^Suematsu e t al.l (|l995l) studied this relation 
and found that longer lived mottles are physically longer and that the relation is 
consistent with the ballistic hypothesis. They noted, however, that although some 
mottles are represented by ballistic motions, a pure ballistic trajectory is unreal- 
istic, because the observed values of initial velocity and lifetime are too small to 
fit the data, while the small dece lerations would requ i re the mottles to be tilted 
typically 60° - 70° from the LOS. IChristopoulou et"all (|200ll) . from light curves of 



mottles, studied their proper motions performing a least-square fit of the ballistic 
motion to the height measurements as a function of time. Their best fit was for 
an inclina tion angle from the local vertical o f ~ 53° and an initial velocity of ~ 
56kms~^. lRouppe van der Voort et al.l ( 20071) used high spatial and temporal res- 



olution observations of a quiet Sun region obtained with SST on La Palma. They 
reported that it was very difficult to track individual dark mottles during their 
lifetime for several reasons among which: a) many mottles lack a sharp top end, 
b) they display significant fading which appears to be induced by opacity changes 
with time and c) many mottles undergo not only up and down motions, but also 
significant motion transverse to the magnetic field. These factors contribute to LOS 
superposition and render the identification and follow-up of mottles throughout 
their lifetime quite problematic. They were able, however, to follow some mottles 
and found that their tops undergo a parabolic path. In addition, they found that 
the deceleration and maximum velocity are linearly correlated and they suggested 
that leakage of global oscillations from the photosphere (with dominant periods 
around 5 min) plays an important role in the formation and dynamics of at least a 
subset of quiet Sun mottles. This finding further suggests that the driving mecha- 
nism for at least some of the mottles is similar to the driving mechanism of DFs. 

Another technique to infer velocities of structures seen in absorption on the 
solar disk has been introduced by Beckers (1964). It is based on the contrast 
profile and is called the "cloud" model (see Section [8. ip . It has been used exten- 
sively in the analy sis of observations of identifiable fine struc t ures in absorption 



on the solar disk (iGro ssmann-Doerth and von Uexkiil l |l97lL Il973l Il977l: iBrav 



1973; Lo ughheadl 11 973: Bray, 1974; Brav and Loug hhead. 1983; Tsir opoula et aL. . 
1994c iTziotziou et al., ,2003) . iGrossmann-Doerth and von Uexkiill (.197l|) . as well 



Bravl(|l973l) . applied this model to chromospheric mottles observed at disk cen- 



tre and derived values for the 4 parameters from the observed contrast profiles. 
iLoughhead ( 19733 , howe ver, showed that this model failed to explain contrast pro- 
files near the lim b, also Brav (Il974l) could not reproduce some contrast profiles 



of sunspot fibrils. l Alissandrakis et al.l (jl99 d) have described the range of circum- 
stances under which this model is applicable and were ab le to map the 2D s patia l 
variation of the obtained parameters. In subsequent work. pTsiropoula et al. I (|1994|) 



using the "cloud model" showed that the predominant pattern of bulk motion in 
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dark mottles is that of down-flows at their foot-points and a lternatinR up-flows and 
down-flows at their tops. This result was re-conflrmed by iTziotziou et alJ (j2003|) 
who used another set of Ha observations obtained by a different instrument. In 
these works the bi-directional, as well as the recurrent character of the velocities 
was revealed. This bi-directional flow followed by down-flowing material along the 
whole structure is shown in Fig. [151 As can be revealed from this figure, the upper 
parts of the mottles show flows in both directions, while the bottom parts show 
flows in the downward direction only and the whole process is repeated with a 
period of ~ 5 minutes. In order to explain this kind of flow pattern they suggested 
that magnetic reconnection might be the mechanism responsible for the formation 
and dynamics of these structures. This suggestion is reasonable, since this process 
seems to be the most suitable physical mechanism that can explain the observed 
flow of material towards and away from the solar surface. Furthermore, this pro- 
cess can occur at the network boundaries where neighboring pairs of flux tubes 
of opposite polarity are driven together by external plasma flows and reconnect. 
Based on these observational finding s, they proposed a simple reconnection model. 
Such a model was also proposed by IPikel'neii (jl969l) . A schematic representation 
of this model is given in Fig. [16] (of course, unipola r as well as bip olar fields of 
appropriate polarity could be involved, see Figure 8 of lWilhelml (|200Cll) ). According 



to this model, the squeezing of opposite polarity field lines at the chromospheric 
level leads to cooling by radiation of the compressed gas, which is trapped between 
them and consequently to a downflow due to gravity (left part of Fig. I16p. When 
opposite held lines come close, reconnection occurs, then part of the material is 
carried upward by the reconnected held lines, while the material below the recon- 
nection region moves downwards under the action of both gravity and magnetic 
forces (middle and right part of Fig. [16}. Th is cycle is repea t ed un til the held is 
annihilated. Analyzing the same set of data. pTziotziou et al.l ( 2004) showed that 



sometimes there is a temperature excess at the location where reconnection seems 
to be occurring, indicating the presence of local heating (Fig. [TSl right). Their 
study of the intensity and Doppler velocity time sequence indicated also that in- 
dividual mottles appear in bursts, lasting ty pically for about 5 min and u sually 
reappear several times at the same location. iMurawski and Zaaarashvilil ( 2010(l 



developed a 2D rebound shock model in which a velocity pulse is steepened into 
shocks which produce mottles/spicules and reported quasi-periodic behaviour of 
the rising material due to consecutive shocks, as well as bi-directional flows due 
to the superposition of falling off and rising plasma portions. 



IContarino et al ] (|2009|l analyzed data acquired along the Ca ii 8542 A line with 
the IBIS instrument at the DST. They applied the cloud model to four mottles 
and obtained values for the four parameters given by this model. They found 
an almost symmetric velocity distribution between —10 and -|-10kms~^, which, 
they suggested, is an indication of the presence of both upward and downward 
velocities inside the mottles. As they reported, however, the LOS velocities they 
found did not show a regular behavior. In one of the mottles examined, they 
obtained alternating upward and downward motion occurring in phase in all the 
examined segments. In another mottle, they obtained up-flows in one endpoint 
and down-flows in the other which reverse with a ~ 4 min period. Finally, in the 
other two mottles, they found both positive and negative velocities which change 
with an apparent irregularity. 
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Fig. 16 A si mple reconnectio n mod el (see text) explaining the observed velocity behaviour 
in mottles, see iTziotziou et alj fcOOaP 



Oscillations in dark mottles Different physical but mainly magnetic conditions in 
the internetwork (IN) and network boundaries (NB) are clearly reflected in the 
oscillatory behaviour of both regions, as well as in dark mottles emanating from 
the network boundaries. Several issues concerning the dyn amics of netwo rk and 
internetwork have been already addressed in the review bv lRuttenI (Il999l) . while 
intensity and velocity oscillations at both IN and NB have been extensively inves- 
tigated by several other authors. 

Nowadays, the general belief is that NB regions show no power of particular 
significance below 5min, in contrast to IN regions that display a dominant period 



of 3min, especiahy in Doppler shifts (IPame et a .1 1984 :'Deubner and Fleckl 



Lites et al.Lll993l:lvon Uexktill and Kneeil . ll995f 



1990; 



Cauzzi et al.. 2000.:,Kriiger et al 



200l | ). However, dif f erent results have also been reported (see the Introduction 



of iTsiropoula et al.l ( 20091) ). stemming mostly from the use of different meth- 
ods or from observati ons i n different lines formed at different atmospheric layers. 
iHansteen et al.l (j2000l) and iBaneriee et al ] (|200ll) for example, have reported net- 
work oscillations probably attribute d to waves, which are produced in bursts with 
coherence time of about 10 - 20 min. lJess et al. 1 (|2009f ) used observations obtained 
by SST investigated a group of network bright points. They reported full-width 
half-maximum (FWHM) oscillations with periodicities ranging from 126 to 700 s 
originating above a bright point. They also reported a lack of cospatial inten- 
sity oscillations and transversal displacements which, as they stated, rules out the 
presence of magnetic-acoustic waves. They suggested that the FWHM oscillations 
could be produced by a torsional Alfven perturbation. 

The li nkage between NB osc illations and mottles has been extensively ad- 
dressed bv ITziotziou et all ( 20041) in a study concerning periodicities of fine struc- 
tures in a quiet Sun network region where a chain of mottles was observed. Con- 
cerning quiet Sun mottles, although differences do exist in the periodicities of 
intensity and velocity variations, which are often bursty and intermittent, with 
velocity variations resulting in a large range of periods from 100 to 500 s, individ- 
ual mottles or regions of mottles do exhibit a most prominent period in the 5 min 
range (see FigjlTp. They also show a 3 min signature which, however, is never the 
domin ant one. Furthermore, a n intermittent signature of 100 s (lOmHz), reported 
also by IHansteen et al. ] (|2000l) . seems to exist; it is, however, unclear whether this 
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short period could be associated with the dynamics of recently reported Type II 
spicules. 



More recently. iTsiropouIa et aL in a multi- wavelength analysis of a solar 



network region showed that both mottles and NB show a periodicity of ^ 5 min in 
all considered lines. Phase differences in the network boundary region indicated an 
upward propagation of waves, while in the region of mottles the phase difference 
was mainly negative for periods of 250 - 400 s suggesting a downward propagation 
due probably to refraction of waves from the inclined magnetic field of mottles. 
A dominant peak around 1000 s has also been reported which is, however, unclear 
if it is related to the presumed in literature lifetime of mottles. The relationship 
between the oscillatory properties of a network region and th e magnetic configu- 
ration of the chromosphere has bee n the sub.ject of stud y by iKont ogiannis et al.l 
([2010b). Based on a previous paper ( Kontogiannis et al.Ll20 10a). they showed the 



existence of power enhancement and suppression (power halos and magnetic shad- 
ows respectively) in 3, 5 and 7 min Hq oscillations and they concluded that p-mode 
leakage, mode conversion, as well as reflection and refraction of waves on the mag- 
netic canopy (see Section [6] for further details) can play an important role to the 
observed properties of network oscillations. 

The dominant 5-min periodicity present in both NB and dark mottles is com- 
patible with the idea of the leakage of photospheric 5-min acoustic oscillations, 
which are the well-known p-modes, to higher layers. Commonly, it was accepted 
that only waves with periods above the acoustic cut-off period (which at the 
chromospheric level is of the order of 3 min) could travel freely in the solar atmo- 
sphere, whereas waves with higher periods were considered as evanescent. However , 
sever al theoretical studies ( Michalitsanosl Il973l : iBel and Lerovl . Il977l : ISuematsui 



1990f ) have shown that the acoustic cut-off period is increased when waves are 
propagating along inclined flux tubes allowing the leakage of p-modes to higher 
l ayers . For the first time p-mode leakage has been suggested bv lDe Pontieu et al.l 
(|2004h to be related to the formation of active region fibrils, p-mode leakage, how- 



ever, seems to be prese nt also in quiet Sun mottles which are also inclin ed fiux 
tubes (Suematsu, 1990i : iTsiropoula et al.l , 120091 : IKontogiannis et all l2010a. b.) . It 
should also be stressed that measured periods of the order of 5 min in IN regions 
close to NB could well be attributed to the presence of inclined mottles that cover 
these regions. Recent results reported by Jefferies et al. (2006) show that a siz- 
able fraction of the photospheric acoustic power at periods above the acoustic 
cut-off might propagate to higher layers within and around the magnetic NBPs. 
They referred to these regions as "magneto-acoustic portals" and showed that 
through these regio ns a significant p o rtion o f energy is provided for heating the 
solar chromosphere. [Khomenko et al.l ( 2008a|) has also suggested that propagation 



of longer-period waves (i.e. equal or greater than 5 min) is even possible in vertical 
magnetic fiux tubes, resulting from radiative losses in NBPs that may lead to the 
lowering of the acoustic cut-off frequency. Recently, iMurawski and Zaaarashvilil 
(|2010[ ) developed a 2D rebound shock model which predicts quasi-periodic raising 



of chromospheric jet-like structures with a period of ~ 5 min. In their simulations 
this periodicity results from a non-linear wake that is formed behind a leading 
pulse rather than from p-mode leakage. They found that the periodicity strongly 
depends on the amplitude of the initial pulse and can be longer for stronger pulses. 
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Fig. 17 Wavelet analysis results for a) the velocity and b) the Ho line-center intensity vari- 
ations along the axis of a mottle, showing a dominant period close to 5 min. Top panels show 
the global power spectrum as a function of period and position along the length of the mottle 
while bottom panels show the variations of the global probability (solid line) obtained with 
a randomization method and the peri od (dotted line) of th e maximum global power peak 
corresponding to this probability (from iTziotziou et al.l | |2004D ) 



4.2.2 Fibnls 

The term fibrils refers to mottle-like, dark (absorbing), highly dynamic, jet-hke, 
on-disc structures which are observed above or in the direct vicinity of active region 
plages. The term fibrils refers also to dark structures which are observed in the 
penumbra of sunspots (penumbral fibrils) radially expanding around the umbra 
and creating a filamentary structure and to structures in active region plages which 
do not show jet-like behaviour; both are probably associated with low-lying loop- 
like structures that connect regions with opposite polarity magnetic fiux. Jet-like 
fibrils related to active region plages are relatively shorter than penumbral fibrils 
or non jet-like fibrils or mottles and probably constitute also a fraction of the 
observed on-limb spicules. 

A close inspection of Ha filtergrams of plages reveals several marked resem- 
blances to the quiet Sun regions. In particular, the basic components of the rosette 
structure, i.e. NBPs and dark mottles are similar in several respects to the fine 
structure of plages, i.e. facular granules and dark fibrils. The principal differences 
of the two regions are: a) the number of facular granules in a plage is much larger 
than the number of NBPs in the centre of a rosette and b) the fibrils around 
a plage are g enerally much more elongated than the mottles around a rosette. 



iFoukall ( 197lh proposed a continuous relationship between the basic constituents 



of the quiet and active chromosphere. According to this author there is a grad- 
ual progression between the mottles oriented more close to the vertical and the 
fibrils which are rather horizontally defiected becau se of the stro nger field which 
constrains the plasma to follow low- lying field lines. Ipbukall ( 197ll) gives some ob- 



served parameters of fibrils, e.g. length, lifetime, density, temperature, velocity and 
magnetic field and compared them to the corresponding parameters of spicules. 
He concluded that the corresponding parameters are rather similar except from 
the length (fibrils are longer) and the magnetic field strength which is higher in 
fibrils. March (1976) performed a systematic study of the lifetimes, lengths and 
evolution of fibrils in an active region. He reported that all fibrils he examined 
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Fig. 18 Ha wideband {left panel) and Ha line center {right panel) images of NOAA AR 10813 
taken by SST. The white and black contours in the left image outline the positive and negative 
magnetic flux, respectively. The DFs are predominantly observed in the mostly unipolar plage 
region between the two s unspots. The arrow indicates the direction of the disk center (from 
iDe Pontieu et al.l ll2007al'l 1 



evolved by an upward extension followed by a downward retraction. He also found 
that the lifetime of a fibril is a monotonically increasing function of its maximum 
apparent length. Based on this relationship, together with the variation of fibril 
lengths with time he suggested that fibrils result from material being impulsively 
injected into magnetic field lines at ^ 30kms~^, and returning back under gravity. 

iPietarila et al.l (j2009) presented high-resolution observations of an active region 
in the Can K line obtained from the SST on La Palma. They found that very thin, 
bright fibrils are a prevailing feature over large portions of the observed field. Fibrils 
have not been observed with such detail in this line before for several reasons. One 
of them is that observations in this line are, usually, obtained with broad-band 
filters. However, only the Can K line core is purely chromospheric, while the wings 
are formed in the photosphere. With filters having a large passband, the intensity 
is summed over a significant portion of the line profile and thus the intensity of the 
line core originating from the chromosphere is rather small relative to the intensity 
of the line wings originating from the photosphere. If, however, observations are 
performed in a sufficiently narrow wavelength band to avoid contamination by 
photospheric radiation then fibrils can be observed. Another reason is that these 
structures are very thin and thus very high-resolution observations are required. 
The fibrils are structured by the underlying magnetic field and the bright endpoints 
are clearly co-spatial with the magnetic concentrations in the photosphere. This is 
not as obvious in the strong plage where the density of the fibrils is high enough 
to make identifying the endpoints difficult. They reported that in a strong plage, 
where there is more magnetic flux and due to unipolar crowding the field is more 
vertical, and so are the fibrils, which are also short, and form a thick carpet 
covering the underlying photosphere. At the plage edges, on the other hand, the 
fibrils become longer and more organized, they extend radially away from the 
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magnetic concentrations, are nearly parallel to one another, more inclined and 
extend over multiple granules towards the quiet Sun. 

While mottles have been the major constituent of the quiet solar chromo- 
sphere, recent high resolution observations revealed that in active regions there 
exist mottle-like structures called dynamic fibrils (DF) which appear to form a 
subset of what have traditionally been called active region fibrils. Another subset 
is constituted from fibrils that do not show jet-like behaviour and which are appar- 
ently low-lying, heavily inclined structures connecting opposite polarity magnetic 
fiux. This subset of much longer fibrils is more stable in time than the DFs. 

Length s, wi dths and lifetimes o f DFs D Fs were extensively studied bv lHansteen et al 



(|2006h and lDe Pontieu et al.l (|2007al) who have measured their morphological prop- 



erties and described their dynamic behavior using very high spatial (120 km) and 
temporal resolution (Is) Ha observations obtained with the SST on La Palma. 
They found DFs are shorter than mottles; their mean length is 1 250 km with 
values ranging between 400-5 200 km. Their average width, which does not vary 
much with time and height, is ~340km with most of them being 120-380 km wide, 
although much wider structures can be found (~ 1 100 km). The lifetime of DFs 
is shorter than that of mottles and lies between 120 s and 650 s, with an average 
of 290 s. These properties are subj ect to regional dif f erences , i.e. they vary signif- 
icantly with position in the FOV. iDe Pontieu et al ] (|2007a[) foun d that DFs (see 



Fig. | 1 8l) were shorter in length and lifetime in the denser plage. IPietarila et al.l 
(|2009() applied a method to isolate individual fibrils and measured their widths 
and lengths. They found that most fibrils have widths between 80 km and 100 km, 
while widths of fibrils in a strong plage have a wide distribution with a max over 
150 km. I'Anan et al l (2OI0I) using SOT/Hinode Can H observations of an active 
plage identified 169 fibrils. They found that 80% were "parabolic" showing a cy- 
cle of rise and retreat and 10% "faded out", i.e. did not show a complete retreat 
cycle. It should be pointed out that these authors use the term "parabolic" (or 
"ballistic" ) motion in the context of constant deceleration during the downward 
phase which is not necessarily a free-fall dominated by gravity. The mean heights 
of the "parabolic" and the "fade out" fibrils were 1 300 km and 820 km, while their 
lifetimes were 179 s and 197 s respectively. It has not been examined yet which is 
the exact relationship between the DFs observed in Ha with the SST and those 
observed in Can H with the SOT/Hinode broadband filter which for on-disk ob- 
servations gives a predominantly photospheric signal. 



Evolution and velocities of DFs iDe Pontieu et aP ( 2007ah studied the evolution of 



these structures using very high spatial and temporal resolution Ha observations. 
In order to find the temporal behavior of individual DFs they constructed x — t 
plots showing the evolution of the extent of a DF as a function of time. They 
analyzed 257 DFs and concluded that the length of the DFs as a function of time 
is well described by a parabola (see, e.g. Fig. ll9p . Fitting the time evolution of the 
length of the DFs with a parabola they derived the maximum velocity (on ascent 
or descent), maximum extent, deceleration and duration. The parabolic paths are 
characterized by a large initial velocity, usually of order of 15 - 20kms^^, that de- 
crease linearly with time. The deceleration is, usually, between 120 and 280 ms^^. 
The deceleration is typically only a fraction of the solar gravity and incompatible 
with a ballistic path at solar gravity. For the particular DF shown in Fig. I19lthev 
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Fig. 19 An X — t plot of the evolution of a DF (along the projected path). It follows a near 
perfect parabolic path in its rise and descent. The white dashed line indicates the best fit used to 
derive deceleration, maximum velocity, duration and maximum length (from lDe Pontieu et al.L 
l2007a^ 



found a deceleration of 216ms^^, a maximum velocity of 27kms^^, a maximum 
extent of 1 800 km, and a lifetime of 240 s. They also found that DFs exhibit lower 
velocities and higher deceleration in the denser plage and that these differences 
cannot be attributed to projection effects. They also constructed scatter plots of 
various DF properties found in two different regions within the FOV. One of the 
regions contained a dense plage, i.e. strong magnetic field concentrations in which 
the magnetic field is generally more vertical, while the other was located at the 
edge of the plage region containing more inclined magnetic fields. They reported 
significant differences of fibril properties between those occurring above the dense 
plage and those occurring at the edge of the plage. From the scatter plots they 
derived some intriguing correlations. One interesting correlation is that between 
the deceleration and the maximum velocity of the DFs which shows a clear linear 
relationship (Fig. 1201 top left). Another interesting correlation is that between the 
maximum length and the duration of the DFs (Fig. [20l top right). The longer DFs 
tend to have longer lifetimes. They also found that the maximum velocity and 
maximum length of DFs are well correlated: DFs with higher maximum velocity 
tend to be longer (Fig. 1201 bottom right). The deceleration of DFs shows a some- 
what less clear correlation with the DF duration (Fig. (20] bottom left). Although 
it seems that longer lived DFs typically experience lower decelerations, there is 
a large spread in the values. There are cle ar differences betw een the two regions, 
while the correlations are not quite linear. lAnan et all (|201Q) found that the DFs 
they examined follow a ballistic motion with constant deceleration. They found 
mean maximum velocities of the parabolic and fade out fibrils of 34kms^^ and 
16kms^^, respectively, and mean accelerations of — 510ms~^ and — 130ms~^, 
respectively. They also found that the deceleration was proportional to the max- 
imum velocity, i.e. the initial velocity of the ejection. The range of the maximum 
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Fig. 20 Scatter plots of deceleration vs maximum velocity {top left), maximum length vs 
duration {top right), deceleration vs duration {bottom left), and maximum velocity vs maximum 
length (bottom right). DFs in two different regions were considered: one region contained dense 
plage {red asterisks) and the other was found at the edge of the plage {blu e diamonds). The 
scatter plots for the two regions often seem to have different slopes (from iDe Pontieu et all 
l2007eS) ) 



velocity was between 6kms - 110 kms for the parabolic fibrils, whereas it was 
between 5kms^^ - 40kms^^ for the fade out spicules. 

Apart from the above mentioned analyses of fibrils based on filtergrams, only 
a few q uantitative spectroscopic studies have been undertaken. iLangangen et all 
(|2008ah obtained spectrograms in the Ca ii 8662 A line using the SST spectrograph. 
They identified 26 DFs and infer Doppler shifts from the spectra. They also de- 
rived mean values ^ 89ms~'^ for the deceleration, ll.Skms^^ for the maximum 
vel ocity and 217 s fo r the li fetimes. These values are lower than those obtained 
by [De Pontieu et al. ('2007a') who measured proper motions in narrow-band im- 
ages. |L;mg£mgenjet_a^ (200 8a) performed also numerical simulations and explained 
the discrepancies in the lower maximal velocities derived from Doppler measure- 
ments compared to the proper-motion velocities as due to both the low formation 
height and the extensive width of the contribution function of the Can 8662 A 
line. They also concluded that their observations support the result that DFs 
are driven by magnetoacoustic shocks excited by convective flows and p-modes. 
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Fig. 21 Wavelet power spectra for Ha-350 mA a nd Ha+350 mA chromospheric oscillations 
in fibrils showing intermittent 5-min periods ffrom lPe Pontieu and Erdelvil ll2006l )) 



iLangangen et al" I (|2008bl) used a time series of filtergrams obtained in the red and 
blue wings of the Ha line near the solar disk center and measured proper motions 
and Doppler signals in 124 DFs. They obtained mean lifetimes of ~ 258 s. Again 
they found Doppler velocities to be a factor of 2 - 3 smaller than velocities de- 
rived from proper motions in the image plane and the corresponding decelerations 
to have a difference of a factor of 5. They attributed the difference to the radia- 
tive processes involved, since the Doppler velocity originates from a wide range of 
heights in the atmosphere, in contrast to the proper-motion velocity, which is a 
very local quantity, because it is measured from the sharply defined bright tops of 
the DFs. 

A combination of high-resolution observations and advanced numerical simula- 
tions have shown that DFs are most likely driven by magneto-acoustic shocks that 
form w hen photospheri c oscill a tions leak into the chromo sphere along inclined flux 
tubes (|Hansteen et al.l (|2006t) : iDe Pontieu et al.l (|2007al) '). The inclination of the 
magnetic field lowers the acoustic cut-off frequency sufficiently to allow p-modes 
with the dominant low frequency to propa gate along flux tubes (Michali tsanosl 
I1973I: iBel and Lerovl[T977HSuemats"uLll990l) . The results of the simulations span- 
ning from the upper convection zone to the corona, lead to the conclusion that 
DFs are formed by chromospheric shocks driven by global p-modes and convective 
fiows. 



Oscillations in fibrils Although quiet and active regions are considered to be made 
up by the same elementary magnetic fiux tubes, there are important differences 
in the observed magnetic flux densities mainly due to the different filling factors, 
which are manifested also in the respective oscillatory behavior. However, far less 
work exists in the literature for oscillations in plages and fibrils compared to work 
for oscillations in IN, NB and mottles. 

iBhatnagar and Tanaka (1972) ex amining Ha filtergrams have reported 5 min 
oscillations in plages, while iMuglachI (j2003l) using UV observations from the Tran- 
sition Region and Coronal Explorer (TRACE) found 5 min oscillations in both 
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Fig. 22 DOT picture of the sun taken on June 18, 2003 in the G band, the Call H wing, and 
the Call H core. The region is close to the solar limb. Scale: the tick marks are separated by 
1". The magnetic network, shows up as bright patches in the Call H wing and the G band. 
The l ong thin struc tures in the Can H core sticking out of the network are termed "straws" 
(from lRuttenI 



plage and network regions and 3min oscillations in IN regions. lOe" Pontieu et al.l 



(|2003h examining high-resolution SST observations of dynamic fibrils in a plage 



region reported the presence of significant oscillatory power, although only a small 
fraction of dynamic fibrils showed oscillations or recurrence. Oscillatory periods in 
dynamic fibrils range between 4 and 6 min, similarly to dark mottles in quiet 
Sun r egions, and seem to be concentrated mostly in the m ore inclined struc- 
tures (jPe Pontieu et al.l |2004 iDe Pontieu and Erdelvl |2006|) . This similarity in 



the oscillatory b ehavi or of fibrils and dark mottles has been clearly shown by 
iTziotziou eraP (|2004l) . who analyzed a plage region containing almost vertical 



fibrils, named as dark grains, which exhibited a unique, high-probability period 
of ~ 5 min. This period has been explained as the period signature of either flows 
along vertical structures anchored in the photosphere or as a possible signature of 
a maximum amplitude p-mode interference wavetrain, compatible with the idea 
of p-mode leakage reported bv IPe Pontieu et al. (2004 ^ or as due to consecutive 
shocks as pre dicted by the 2D r eboun d shock of Muraw ski and Zaaarashvili ( 2010). 
Furthermore. pTziotziou et al.l ( 20041) showed that the fibril-free area, although in 



close spatial proximity to the fibril area, had IN-like characteristics with periods 
closer to 3 min in contrast to the 5-min period of fibrils which is similar to the 
period obtained for dark mottles. 



The oscillatory characteristics of DFs, as IPe Pontieu et al.l ( 20041) have shown 



exhibit similarities to photospheric p-mode oscillations, such as variable periods 
between 200 - 600s and intermittency (existence of wave packets). Moreover, ob- 
servations have indicated that DFs tend to occur near the per ipheries of plages 
where magnetic field lines are often more inclined. This has led iDe Pontieu et al.l 
f2004) to suggest that DFs are driven by l e akage of normally evanescent p-modes 
in inclined magnetic fields, as ISuematsul ( 1990|) had also previously suggested. 



something that can explain the observed 5-min periodicities. 
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4.2.3 Straws 



iRuttenI (|2006l l2007l) using high-resolution ground-based observations obtained 



with DOT in the wide-passband (of 1.4 A) Can H filter observed some new fea- 
tures, which are finer than the traditional mottles and were termed "straws" 
fFig. [22)) . These structures are long, highly dynamic, bright against dark inter- 
network background, extend radially outward from network bright points, and 
occur in "hedge rows" . The hedge rows of upright straws are bright and thin in 
Can H, bright, optically thick, twice as high and much thicker than in Can H 
in Ly-a, much less distinct and dark in Ha fine cen t er, an d more prominent, but 
less upright and very dark in the Ha wings. iRuttenI (|2006l) concluded that straws 



reflect transition region conditions and the differences in the different lines are 
due to radiative transfer effects. Such features were also recen tly observed with 
the Sol ar Optical Telescope at the solar limb in Can H images. IPe Pontieu et all 
refers to these straw-like features as "Type H" spicules. 



4.2.4 Rapid blueshifted excursions (RBEs) 

A class of short-lived events observed on the solar disk with IBIS /DST and char- 
acterized by large Doppler shifts that appear only in the blue wing of the Ca 11 IR 
line have been recently reported by Langange n et al.l (HoOSc.) . They were found at 
the edges of the rosettes as sudden broadening of the line profile on the blue side 
of the line without an associated redshift and were denoted as rapid blue-shifted 
excursions (RBEs) by the same authors. For the measured RBEs they derived an 
average length of 1.2 Mm, an average width of 0.5 Mm, a mean lifetime of 45±13s, 
and velocities of the order of 15 - 20kms^^. Using Monte Carlo simulations, they 
showed that the observed blue-shifts can be explained by a wide range of spicule 
orientations combined with a lack of opacity in the upper chromosphere. Similar 
events, but longer lived and spatially more extended, have been reported earlier 
( Wang et al.L fl99^ . These events were also revealed using spectral imaging data 



in the Can 854.2 nm and Ha lines obtained with t he CRisP Imaging Spectropo- 
larimeter at the 1-m aperture SST on La Palma ( Rouppe van der Voort et al.l 



|2009[) . These authors used a manual detection scheme together with an auto- 
mated algorithm to detect RBEs. They identified 413 features in the Can 8542 A 
data set and 608 in the Ha data set and claimed that the number of RBEs if 



they were observed at the limb would be 1.9 per linear arcsec. lJudge et al 



exrapolating the number of RBE detections of lRouppe van der Voort et al 



2011) 



20091) 



estimated that there should be ~ 10"' RBEs on the Sun at any given moment. They 
found rapid blue-ward excursions in the line profiles of both chromospheric lines 
(Fig. I23p and suggested that these structures may be the on-disk counterparts of 
Type II spicules. They measured lengths, Doppler velocities and widths of RBEs 
and found that the values of these quantities are consistently higher for the Ha 
RBEs compared to the Can RBEs. Thus Ha RBEs are on average longer (of order 
3Mm vs 2Mm), higher widths (13kms~^ vs 7kms^^) and have higher average 
Doppler shifts toward the blue (35kms^^ vs 15kms^^). 
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Fig. 23 Properties as a function of position of an RBE in the Ca II 8542 A data set (top panels) 
and the Ha data set (bottom panels). Color-coded- measured Doppler velocity (left panels), 
width (middle panels), and mean spectra over three portions of the RBE (right panels: closest 
to foot-point in blue, middle part as dashed green, part furthest away in red and mean spectrum 
in black. Upper curves show the spectral profile, the lower curves show the subtraction of the 
average spectrum and the spectral profile). The extracted RBE is shown with a thin colored 
line in the left and middle panels with blue for the third of the length closest to the foot- 
point, green for the middle part, and red for the part furthest away from the foot-point. The 
measured parameter is shown displaced to the left of the RBE. The background image is at line 
center (middle panels), at a blue position of 36kms~^ (bottom left panel), and a Dopplergram 
at 30kms~^ (top left panel). At both the bottom and top end of the RBEs the amount of 
blue- ward absorption drops significantly. When this happens, the Doppler velocity and width 
is set to zero. This is why the black color coding at both ends of the RBEs is not indicative of 
lower velocities or widths, but rather sign ifies the spatial extent of the blue absorbing feature. 
From lRouppe van der Voort et al.l j2009l'l 



5 Relation between mottles/fibrils and photospheric bright points 

Multi-wavelength photospheric and chromospheric observations of quiet solar re- 
gions, such as the ones presented in Fig. [2l show a clear spatial relationship 
between the dark elongated structures, called mottles, emanating from the net- 
work boundaries and the photospheric bright points. Hence, it is quite natural 
to think that these structures are somehow related. Bright points, which can 
be regarded as proxies of small-scale magnetic elements in the photosphere, are 
prominent at network boundaries in the Hq line-wing images at high spatial reso - 
lution fe.g.'Punn and Zirkeil ll973HSuematsu et allll995tlLeenaarts et al.Ll2006h . 
As[^ematsu et al. (1995) showed, by comparing Ha-0.65 A and Ha-|-0.65 A filter- 
grams, bright points are very prominent at network boundaries i n the blue-wing, 
while only a few appear in the red-wing or line center. Recently, iLeenaarts et al.l 
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()2006[ ) concluded that in the blue Ha wing bright points represent a promising 
proxy magnetometer to locate isolated magnetic elements. 



It is of great interest to investigate whether the appearance of bright points, 
which represent intense magnetic concentrations, is related to fine structure for- 
mation at the network boundaries. However, it i s extremely d ifficul t to directly 
associate bright points to specific dark mottles. IPunn et al.l ( 1974) . used their 
excellent H a filtergrams, but they were not able to find a definite relationship be- 
tween them. lSuematsu et al.l ( 199 5*) . using time series of Ha observations were able 
to connect a bright point with a mottle through its whole life. They reported, how- 
ever, that there were many mottles which were not associated with bright points. 
They also noted that bright points often appear at the foot-points of mottles at 
Ha-0.65A, but during their peak extension or falling phase rather than at their 
beginning. In Ha+0.65A bright points seem to be obscured by falling material. 
In the blue wing, however, and in the falling phase of mottles one can see a much 
deeper layer due to an effective red shift of the line. Enhancement of the brightness 
in —0.65 A could occur due to compression of the atmosphere by in-falling spicule 
material that penetrates deeper the solar atmosphere, delivering its potential en- 
ergy. The infalling material may even distort and kink the fiux tube inward; this is 
inferred from the inward offset of the spic ule base in the r e d-wing com pared to its 
blue-wing position (see Suematsu, 1998: Para et al.l . ll998f ). [Cauzzi et a l. (2009), in 
simultaneous Ha and Ca ii infrared observations obtained from IBIS noticed that 
in Ha at minimum intensity dark structures show sometimes marked brightenings 
towards one end. Comparison with magnetograms showed that these bright end- 
ings correspond to the stronger magnetic concentrations, although the magnetic 
network, which is easie r to see in the respect ive Can infrared images, is not easily 
recognized in this line. ISterling et aD (|2010l) used Can H observations of a north 
pole coronal hole obtained with SOT/Hinode and investigated the roots of some 
spicule-like features seen on the disk near the limb. They found that these fea- 
tures seem to emanate from fast moving Can brightenings. Frequently the ejected 
spicules, which have speeds ~ 100 km s^^, seem to occur when these brightenings, 
which have horizontal velocities of a few lOkms"^ appear to collide and disap- 
pear. These spicules could be associated to the so-called Type II spicules. Based 
on their findings, they suggested that Can brightenings could represent acoustic 
shocks or fast-moving interacting magnetic elements. As they stated it is hard to 
understand how the energy from acoustic shocks could be converted to spicules of 
such high velocities. If, on the other hand. Can brightenings represent different 
magnetic elements, their merging could lead to magnetic cancellation and recon- 
nection, resulting to the deposition of thermal energy in the low atmosphere. 



The existence of a relationship between bright points and the correspon ding 
struct ures found in active regions, called fibrils, has been examined by P ietarila et al.l 
(|2009l) in Can K line observations. They reported that bright foot-points of fib- 
rils are clearly co-spatial with the photospheric magnetic concentrations. As they 
noticed, however, in a strong plage it is not easy to identify such a correlation, 
because the density of the fibrils is very high, the fibrils are more vertical and form 
a carpet covering the surface. 
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6 Mottles and fibrils and the formation of the magnetic canopy 

The highly ordered photospheric magnetic fields of the plages and magnetic net- 
work spread out and occupy the whole atmosphere as they extend upwards. Thus, 
instead of the bright points that constitute the network and the facular gran- 
ules that constitute the plages, which more or less represent the cross-sections of 
rather vertical flux tubes, thin, elongated, and usually inclined structur es domi- 
nate i n the chromosphere and provide evidence of its complex structure. iGabriell 
(|1976[ ) proposed that the magnetic field of the upper solar atmosphere has a wine 



glass geometry as a result of magnetic field lines fanning out of the network. Such 
a general config uration of the chromosphe ric magnetic field was already assumed 
before (see e.g. IKopp and Kuperud . Il968l) and is physically anticipated. Due to 
the almost exponential decrease of the gas pressure with height, the magnetic fiux 
tubes expand and form a funnel-like magnetic geometry with almost horizontal 
fields. One interesting quantity describing the relative importance of gas and mag- 
netic field is the plasma /3 parameter (where j3 = ^.irp/B^). The /? = 1 layer, i.e. 
the layer where the gas and magnetic pressures are equal (or, equivalently, the 
sound speed equals the Alfven speed) is called magnetic canopy. This layer is of 
critical importance because it partitions the atmosphere into contiguous volumes 
of gas pressure (high-/?) and magnetic pressure (low-/3) dominated plasma. The 
location of the magnetic canopy depends on the photospheric magnetic flux and 
on the pressure profile of the atmosphere in which the flux tubes are embedded. 
Therefore, in principle, one may expect that it should be lower near active re- 
gions and higher around the chromospheric network in quiet Sun. In the quiet 
Sun the canopy is situated somewhere between ~ 1 and 2 Mm above rsooo = 1 



(jSolanki and Steinerl . [l990l) . Thus it is within the chromosphere that /3 falls below 



1 and magnetic forces start to gain control on the dynamics of the solar plasma. 
This results not only in mode conversion, refraction and reflection of waves, but 
also, together with the local thermodynamic properties, in the plethora of fine 
structures that characterize this part of the atmosphere. 

Concerning wave propagation, it has been shown by numerical simulations 
and theoretical studies that the canopy acts as a boundary. Waves that propagate 



tion and reflection when they reach the canopy (see e.g. 


Rosenthal et al.l I2OO2I: 


Boffdan et al.'.'2003:'Schunker and Callvll2006l:ICallvll2007l: 


Khomenko et al.ll2008bl 


Kuridze et al., 2009: Nutto et al., 20121). Through the interaction between the var- 



ious MHD modes at the canopy, energy is transferred at higher atmospheric layers. 
This interaction has also been revealed by observational studies which show the 
distribution of the oscillatory power on the FOV. It was found that high frequency 
oscillatory power at the photosphere is increased around inten se magnetic concen- 
trations, in active regions, forming the so-called power halos (iBraun et al. , 1992t 



Brown et allll992t Hindman and Brownl . [l998l:lThomas and Stanchfieldl.l2000l:ljain and Habeil. 
2OO2I: iMudachl '2003V It was also noted that the po wer is reduced a bove strong 



magnetic concentrations at the chromospheric level. iMuglach et al.l ( 2005|) com- 



bined power maps and magnetic field extrapolation in an active region and showed 
that closed lines are associated with a further increase in power, due to reflec- 
tion of acoustic waves. The same areas of enhanced power were less extended 
around nearly vertical open field lines. It is, generally, accepted that this en- 
hancement is due to the interaction of acoustic waves with the canopy. On the 
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Fig. 24 Scatter plots of the average power of Doppler signal variations in a logarithmic scale 
vs the Ha±0.35A minimum intensity {left column) and of t he HaifcO.TA minimum int ensity 
{right column) in a rosette region consisting of mottles from iKontogiannis et ahl l l2010al )) 



other han d, over the n etwork regions the presence of reduced power has been re- 
vealed bv lJudge et al.l (^001) through TRACE and SUMER observations. These 
authors used the term magnetic shadow to describe this power deficit observed in 
areas adjacent to the NBPs and not only above them. Ma gnetic shadow s were 
also detected in chromosph eric lines, such as Can IR or Ha ( Vecchio et al.L [20071 : 
IKontogiannis et al.l l2010a|) . It was explained that the interaction of waves with 
the magnetic canopy is also responsible for the f ormatio n of the magnetic shadows 
(jMcIntosh et al.l l2003h . In IKontogiannis et al.l (|20104 it was further found that 
around NBPs and over rosettes magnetic shadows were detected at Hq±0.35A, 
while at HaiO.TA acoustic halos were detected at the same positions. They also 
noticed that, interestingly, the power maps show a fibrilar structure which corre- 
lated very well with the positions of dark mottles. 

The different types of fibrilar structures seen in chromospheric lines (espe- 
cially in the Ha line) chart inclined magnetic flux tubes. These flux tubes contain 
more atoms and ions in the right state for absorption (or emission) in the corre- 
sponding lines than their surroundings. It is reasonable to assume that the chro- 
mospheric structures are somehow connected with the format ion of the magnetic 
canopy and should play an important role in wave propagation. IKontogiannis et al.l 
( 2010ah examined the relation between power enhancement or suppression within 
a rosette region and the positions of dark mottles. They gave scatter plots (shown 
in Fig. I24p between the average power of the oscillations of the Doppler signal 
and the minimum intensity in Ha±0.35 A and Ha±0.7A (which are obviously re- 
lated to dark structures). In the left scatter plot lower power values correspond 
to lower minimum Ha±0.35A intensity values. The inverse is observed for the 
Ha±0.7A intensities (right scatter plot). The authors arrived to the conclusion 
that power enhancement/suppression is directly related to mottles. Potential mag- 
netic field extrapolation using photospheric magnetograms has shown that mottles 
follow, in most cases , the m agnetic field lines of the chromospheric field (see [25]), 
(jKontogiannis et aI.L r2010bf ). 
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Fig. 25 Ha hne centre 30-minute average intensity image. Overlayed are the field hnes of the 
extrapo lated magnetic field using the potential field approximation (from lKontogiannis et al] 
j2010bl ')'l 



In this same study the authors have assumed a VAL 3C model for the so- 
lar atmosphere and used SOT/SP magnetograms. They obtained the LOS and 
transversal components of the magnetic field, its inclination to the vertical, as 
well as the plasma-/3. From these parameters they calculated the height of forma- 
tion of the magnetic canopy. They showed that the magnetic canopy is directly 
related to the dark mottles, which seem to provide the loci where wave mode con- 
version, refraction and reflection occur. They also showed that the modulation of 
the oscillatory power is defined by the inclination of the magnetic field and the 
relative positions of the magnetic canopy and th e height-of-formation of th e line 
used for the observations. In a subsequent paper [Kontogiannis et all (|201ll) used 
high-precision measurements of the magnetic field in a network region provided 
by SOT/SP on-board the Hinode mission and MDI/SoHO measurements together 
with potential magnetic field extrapolation. They showed that the data of the 
latter instrument lead to higher formed and less horizontally extended canopies. 
The authors concluded that the sensitivity of the instrument used to measure the 
photospheric magnetic field is very important in the derivation of the height of 
formation of the magnetic canopy. 
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7 Trends in observation 

The study of the fine and highly dynamic solar structures requires high angular, 
temporal and spectral resolution observations. Solar observations are moving into a 
new era of high-resolution observations. Ground-based observations, traditionally 
restricted by seeing conditions, are now profiting from adaptive optics and post- 
detection image restoration techniques to achieve a much better angular resolution 
than was possible before, while new large-aperture telescope projects have already 
been initiated. Space-based observations are providing large volumes of data in the 
visible, ultraviolet and X-ray wavelengths, covering extended temporal durations, 
and have revolutionized solar science during the past years with successful missions 
such as YOHKOH, SoHO, TRACE, STEREO, HINODE and more recently SDO, 
while new exciting projects are under way. 

Ground-based solar telescopes are equipped with state-of-the art instruments 
and CCD cameras. The Telescope Heliographique pour I'Etude du Magnetisme 
et des Instabilites Solaires (THEMIS), the 90-cm aperture French-Italian solar 
physics facility in Tenerife, is equipped with a Fabry-Perot filter magnetograph 
(MTR) and a large double-pass spectrograph (MSDP) for multi-line spectropo- 
larimetry and 2-D spectrometry. The 45-cm Dutch Open Telescope (DOT) pro- 
vides solar images at nearly 0.2" resolution while the SST has been recently operat- 
ing a new CRisp Imaging SpectroPolarimeter (CRISP). The Interferometric Bldi- 
mensional Spectrometer (IBIS), a bi-dimensional spectrometry instrument based 
on a dual Fabry-Perot interferometric system, plus ROSA (Rapid Oscillations in 
the Solar Atmosphere camera) are already operating at Dunn Solar Telescope 
of the Sacramento Peak Observatory. The German VTT telescope operates with 
spectropolarimeter instruments such as POLIS (for the iron line pair at 630.2 nm) 
and TIP (for the near infrared) and TESOS which is a 2D spectrometer. The 1.6- 
meter New Solar Telescope (NST) at Big Bear Solar Observatory (BBSO), uses a 
large set of new instruments such as broad-band filter imagers (e.g. a TiO filter 
centered at a wavelength of 705.7 nm, with a bandpass of 1 nm) to deliver high res- 
olution images of small-scale photospheric features with high cadence exposures. 
Knowledge acquired from past and present instrumentation will be used in new so- 
lar telescopes like the 1.5-meter German Gregory Coude Telescope in Tenerife and 
future large-telescopes like the US Advanced Technology Solar Telescope (ATST) 
and the European Solar Telescope (EST). Post-detection image restoration, al- 
though it requires a large amount of processing, enhances the spatial resolution 
close to the theoretical telescope diffr action limit through t echniques, such as the 
phase-diverse image registration (e.g.,'Loefdahl et al.l, fl99i) which is used in SST 
and speckle image reconstruction (e.g., Siitterlin et al., 200lj). An extension of the 
first image restoratio n technique called M ulti-Object Multi-Frame Blind Decon- 
volution fMOMFBD. Ivan Noort et al.l , [20051) operates now at SST. Speckle image 
reconstruction works by dividing the observed field in isoplanatic subfields which 
are reconstructed independently by describing a seeing-sampling sequence of per- 
turbed images called "speckle burst" with a model of atmospheric turbulence, and 
finally rejoining them in one single image. This technique has been extensively 
used for post-restoring of images from several telescopes, such as the DOT, VTT, 
DST and SST. However, although these techniques produce high resolution solar 
images, usage of such images for the quantitative determination of different physi- 
cal parameters with inversion techniques (e.g. cloud modelling) can sometimes lead 
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Fig. 2 6 Left: Deconvolved Hinode SOT image, Right: Reconstructed DST image llWoeer et al.l 



to spurious results and hence possible effects and artifacts of image restoration on 
the data should always be examined before such a scientific analysis. 

Real-time image correction in ground-observations is nowadays achieved with 
adaptive optics. Adaptive optics try to correct the incoming wavefront, similarly 
to the speckle reconstruction principles, by dividing the aperture plane into sub- 
apertures (ideally smaller than the seeing) and applying a wavefront correction to 
compensate for deformations caused by seeing conditions and imperfections in the 
optical system. Adaptive optics makes it feasible to build very large meter(s)-class 
innovative solar telescopes such as the NST at Big Bear Solar Observatory and the 
GREGOR in Tenerife and the even larger future telescopes like ATST and EST. 

Adaptive optics systems are now commonly used on several advanced solar 
telescopes to enhance the spatial resolution of the recorded data. Post image re- 
construction techniques such as speckle interferometry can be used to achieve 
near-diffraction limited resolution over a large FOV rather than the limited area 
where the wavefront sensor of t he adaptive optics system measures the wavefront 
distortions. IWoeer et al.l ( 2008j) investigated the reconstruction properties of the 
Kiepenheuer-Institut Speckle Interferometry Package (KISIP) code. Comparisons 
were made between simultaneous observations of the Sun using the ground-based 
Dunn Solar Telescope and the space-based SOT/Hinode telescope. The analysis 
showed that near real-time image reconstruction with high photometric accuracy 
of ground-based solar observations is possible, even for observations in which an 
adaptive optics system was utilized to obtain the speckle data (Fig. I26[) . 



8 Trends in interpretation of line profiles 

The on-disk solar chromosphere can be observed in the cores of strong absorption 
lines at visible or near-UV wavelengths, while many of the emission lines in the 
EUV part of the spectrum arise in the upper chromosphere. Fine-scale chromo- 
spheric features are usually observed in the Ha Lya, Can H & K and the Can 
infrared triplet. In principle the solar chromosphere is optically thick at these wave- 
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lengths and, consequently, their computation from a model atmosphere requires 
solving the well known radiative transfer equation together with the statistical 
equilibrium equations self-consistently, a problem which is both highly non-linear 
and non-local. However, a simple method called the cloud model, which does not 
involve such elaborate computations and assumes optically thin structures and a 
constant source function along the observed structure has been successfully applied 
over the y ears for several chromospheric structures observed in Ha (see review by 
iTziotzioul i2007) for the determination of physical parameters. Below we discuss 
in more detail the cloud model, other radiative transfer techniques and recent and 
future trends in interpretation of line profiles. 



8.1 The cloud model 

For the deduction of different physical parameters of chromospheric structures 
from observations in the Ha l ine a method i s used whi ch is based in the so calle d 
"cloud" model, introduced bv lBecker^ ( 19641) . (see also lAlissandrakis et al. l ll990D . 



The model is valid for optically thin structures located well above the background. 
It considers for an observed intensity profile I{AX) the contrast profile: 

with a Gaussian wavelength dependence for the optical thickness: 

r(Z\A) = roe"^ ^^d > (5) 

where Io{A\) is the reference profile emitted by the background and AXi = Xqv/c 
is the Doppler shift with Aq being the line center wavelength and c the speed of 
light. The four adjustable parameters of the model are the source function S, the 
Doppler width AXd, the optical depth tq and the line-of-sight (LOS) velocity v. 
All these parameters are assumed to be constant through the structure. 

Once the four parameters are determined from the cloud model, the estimation 
of several other parameters can be made. From the calculated Doppler width values 
and if we assume a value for the micro-turbulent velocity, ^t, we can deduce the 
temperature, Tg. From the optical depth at line center, which may be written as: 

TO = -X^'^ (6) 

and AXjj, N2, e.g. the number density in the second hydrogen level, can be ob- 
tained. Then the electron density, Ne, the total particle density of hydrogen (i.e., 
neutral plus ionized), N^, the gas pressure, p, the total column mass, m, the mass 
density, p, can be determined from the following relations: 

Ne = 3.2 10® ^/n^ cm~^ (7) 
iVjj = 5 10*10°-^'°s^^ (8) 



p=k{Ne + im5lNH)Te 



(9) 
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m = {NnmH + 0.0851Nh x3.97mH)d (10) 

(11) 

where d is t he path length along the LOS (for details about the derivation of these 
parameters iTsiropoula and Schmieded (|l997l) '). 

8.1.1 Optically thick radiative transfer 

For any given species the radiation field is determined from the source function. 
The source function is derived from the atomic level populations that are in turn 
determined by the local radiation field. This inter-dependence defines the highly 
non-linear and non-local aspect of optically t hick radiative tran sfer. 

As already mentioned (see Section I4.2.1|) iGiovaneilH ( 1967al lb[) did pioneering 



extensive NLTE computations of the Ha line contrast profiles relevant to chro- 
mospheric fi ne-scale modelling assuming ID slabs illuminated by the surrounding 
atmosphere. iHeinzel and Schmiederl ( 19941) used Ha line profile observations of 



bright and dark mottles together with a grid of NLTE models of prominence-like 
structures considered as vertically-standing ID slabs irradiated from both sides 
by an isotropic incident radiation and derived the physical conditions in these 
structures. 

In general, most radiative transfer codes used so far for the diagnostics of iso- 
lated solar structures, combine the accelerated (or approximate) Lam bda iteration 
(ALI ) or multilevel ALI (MALI) schemes (e.g.lCannon . 1973; Scharm er and Carlssonl 
19851: iRvbicki and Hummed. Il99ll Il992l: lOlson et al.l Il986f ) wi th efficien t algo- 
rithms for computing the radiation field (Feau trier's method , |Feautried, Il 964 : 



Aued . fl967h and direct integration techniques fe.g. lMihalas et ali ri978: Olson and Kunas j . 
19871) in order to solve the intricate radiative transfer problem. 



There alre ady exist sever al ID-geo metry radiative co des, such as MULTI ( Carlssoiii 



19861) . MALI (|Heinzeifl995l) and RH (|UitenbroekLl200lL the ID version) that com- 



pute the radiative transfer problem in plane parallel atmospheres, using atoms 
containing several or even hundreds of levels and/or thousands of lines. However, 
the solar chromosphere has a complex 3D structure which should be taken into 
account when evaluating the radiation field with non-Local Thermodynamic Equi- 
librium (NLTE) codes. Such codes are usually required to use local ALI operator, 
since any non-local operator would require computationally expensive inversions 
of very large sparse matrices. Since the typical size of 3D radiative transfer prob- 
lems is quite large, some form of parallelization (either over frequency and an- 
gle, or using domain decomposition, or a not-yet implemented hybrid paralleliza- 
tion strategy) is required to speed up convergence to a physical solution. Nowa- 
day s, there exist three 3D NLTE radiative transfer codes, the MUGA code (see, 
e.g [Fa biani Bendicho et al.L Il997l : iBotnen and Carlsson'. 1999") the RH code by 
Uitenbroek (2001) and the Multi3D~ (|Fabiani Bendicho and Truiillo Buend . Il999l: 



Leenaarts et al.t 2009() code, for calculating chromospheric models. 



The situation that prevails in the atmosphere of the Sun poses several dif- 
ferent, problems in the computation of the radiation field with NLTE codes. For 
example, the hydrogen line radiative transfer problem, requires to take into ac- 
count non-equilibrium hydrogen ionization in the chromosphere and hydrogen level 
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populations out of statistical equilibrium. This stems from the fact that the as- 
sumption of statistical equilibrium for hydrogen fails in the chromosphere, because 
the timescale of hydrogen ionization and recombination is of the same order as the 
hydrodynamic timescale. The existing codes fail to reproduce the Ha line center in- 
tensity, as well as the fibrilar structure of the chromosphere. This is because apart 
from taking into account non-equilibrium ionization processes, 3D NLTE time- 
dependent radiative transfer codes are required in order to compute the detailed 
hydrogen line profiles and get reasonably realistic electron density and temperature 
which are very important parameters for the modelling of the chromospheric struc- 
tures. On the other hand, the Can lines suffer less than Hq from non-equilibrium 
ionization processes and thus can be modeled assuming statistical equilibrium. 
The cores of Can H & K, however, form in partial redistribution, which is compu- 
tationally far less trivial than complete redistribution and still has not been fully 
implemented in 3D. To make the problem of simulating these lines even worse we 
mention that the Can H & K lines are usually observed with wide filters which 
deliver intensity contributions from a wide range of formation heights. Partial re- 
distribution effects are less important for the Can infrared triplet and, moreover, 
these lines can be computed in 3D. Current state-of-the-art instruments have de- 
livered many high-quality chromospheric observations in the Can 8542 A line, in 
particular, although this line has a lower line-core opacity than Hq and Ca ii H & : 
K. The NLTE Can 8542 A c omputation using a 3D simulation ('Leenaar ts et al.i 
I2OO9I : iLeenaarts and Carlsso n. 2009), does not show, however, the fibrilar elon- 
gated structures delineating magnetic field lines around the network. This lack 
may be do to several reasons, such as the low resolution of the grid used for the 
simulation, need for larger scale magnetic field configuration, etc. Another pos- 
sibility is that Can 8542 A opacity in the fibrils is caused by physical processes 
curre ntly included in the exi sting simulation codes. On the other hand, quite re- 
cent lylLeenaarts^eT^l] ( 2012f ) using 3D NLTE radiative transfer computations and 



MHD simulations investigated the Hq line formation in the solar chromosphere. 
They were able to reproduce the fibril-like dark structures seen in Hq line-core 
observations. Furthermore, their simulations support the commonly held notion 
that fibril-like structures indeed trace the magnetic field. 

Knowledge of the physical conditions in isolated structures is very important in 
order to elucidate the physical processes which are occurring in them. A state-of- 
the-art method to improve this knowledge is the forward modelling by which data 
from a combination of radiation-MHD simulation codes with 3D NLTE radiative 
transfer codes are converted into observed quantities and after being compared 
with observations can give more reliable initial values for the codes. Although 
existing codes still suffer from severe limitations, they are producing promising 
results, e.g. the computation of the emerging Can 8542 A and Hq line profile s 
with a 3D NLTE model using an RMHD simulation (jLeenaarts et al.l . l2009l.l2012l) . 



the spe ctral line response to r econnection events in the photosphere and transition 
region ( Heggland et al.l . l2009l) . Monte-Carlo simulations of the line profile response 



to upward propagating m aterial, such as in the case of rapid blue-shifted events, 
(jLangangen et an . l2008ch . and toy models of the solar atmosphere to qualitatively 
explain the appearance of the solar limb in the Hinode Ca n H passband caused 
by la rge line-of-sight velocities or wide spicule line profiles ( Judge and Carlssonl 
I2OIOI) . 
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9 Inter-relations between the different fine-scale chromospheric structures 
seen on the disk and at the limb 

One of the outstanding questions in spicule research concerns the relationship be- 
tween spicules at the limb and structures observed on the disk. Although spicules, 
mottles and fibrils are all nowadays treated as jet-like features of the solar chromo- 
sphere and the similarity of many of their properties strongly suggests that some 
of these phenomena are related, if not the same, there has been a strong debate 
in the solar community about the exact relationship between them. In particular, 
it has long been suspected and even assumed that the dark mottles are the direct 
manifestations of spicules, and considerable effort has been expended to prove this 
relationship. The mass velocity is one of the p roperties that could help estab- 
l ishing a kinship between mottles and spicules. iGrossmann-Doerth and Schmidtl 



(jl992i) have examined two distributions of axial velocities, one based on limb the 
other on disk observations. They arrived to the conclusion that the two distribu- 
tions disagree, that the velocity in spicules is significantly greater than in mottles 
and thus that mottles are not the disk counterparts of spicules. They discussed 
several possibilities for this discrepancy, but the most important we think — not 
discussed in that paper - is that the velocities from which these distributions were 
derived have been obtained by different methods: for the limb velocities distribu- 
tion the Doppler shift method was used, while for the disk velocities distribution 
they used the "c loud" model. These two techniques yield very different results. 



as 



it is stated bv lTsiropoula et al.l (|l994l) . The discrepancies on the velocities be- 



tween structures seen on the disk and at the limb, apart from the problem in 
associating with spicules some easily identified class of disk structures, must, fur- 
thermore, be sought as due to several other reasons. These include various selection 
effects arising from different visibilities of the features, different wavelengths used 
for their observations, different phases in their evolution and also eventual bias 
in the data arising from the methods of analysis. They are also arising from the 
fact that the structures are observed under very different conditions: in one case 
(limb) at great heights and along the horizontal LOS; in the other case (disk) at 
low heights and along the vertical. Direct observational evidence for the identifi- 
cation of fine d ark mottles with spicules w ould be that of a dark mottle crossing 
the solar limb. IChristopoulou et al.l ( 200lh found several examples of individual 



mottles crossing the solar limb and provide further support to the association 
between spicules and mottles. On the other hand, some authors suggested that 
similarities in dimensions, lifetime, and physical conditions, such as density and 
temperature, provide indirect evid ence to the identity of spicules and dark mottles 
(jTsiropoula and Schmiedeii [l997l) . 



It is also not clear which the relationship is between quiet Sun dark mottles and 
active region long fibrils and short dynamic fibrils (DFs). Long fibrils are found at the 
edges of plages, are longer, more stable, longer lived and more horizontally oriented 
than the most dynamic DFs, which are found towards the center of the plages, are 
more vertically oriented and have shorter lifetimes (although their lifetime depends 
on their inclination being longer for the more inclined ones). The more complex 
topology of the quiet Sun with its mixed polarity magnetic fields are to be opposed 
to plage regions which are considered as unipolar with magnetic fiux tubes packed 
close together. Due to these differences the conditions between quiet Sun and plage 
regions should differ considerably and this must be refiected in the dynamics, as 
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well as the mechanisms responsible for the formation of mottles and fibrils. Thus 
for mottles magnetic reconnec t ion has been suKgested as th e more likely formation 
mechanism ( Tsiropoula et all Il994l : iTziotziou et all |2003|) , while for fibrils mag- 
netoacoustic shock waves driven by leakage of convective motions and oscillations 
from the photosphere into the chro mosphere along magnetic field s have been sug- 
gested as the formation mechanism. iTsiropoula et alJ ( 19941) . and ITziotziou et al.l 
derived LOS velocities along mottles using the cloud model. Examining the 
temporal evolution of the LOS velocity they found that the predominant motion 
is downward at their foot-points and alternatively upward and downward at their 
tops. Mottles are found at the network boundaries where bipolar elements are 
swept by the supergranular flow. Interactions of the magnetic fields have as a re- 
sult the enhancement of the fiux concentration in the case of same polarities or its 
cancellation in the case of opposite polarities. The bipolar nature of the flux cancel- 
lation process most likely involves magnetic rec onnection and thus it is reasonable 
to consider this process as the driver of mottles. iMurawski and Zaqarashvilil ( 201o|) 
studying the formation of chromospheric jet-like structures in the framework of 
the rebound shock model which was based on a localized velocity pulse launched 
from the lower atmosphere which quickly steepens into shocks. They were able 
to explain the observed speed, width, heights and periods of type I spicules, as 
well as the observed multi-structural and bidirectional flows. Fibrils, on the other 
hand, are found at active region plages which are considered as unipolar, with 
high filling factor of the magnetic field, while the field lines are only moderately 
inclined and the spatial distribution of magnetic flux in their core evolves slowly. 
It is possible that different mechanisms have to be considered as their driver. For 
dynamic flbrils examination of their proper motions have shown that their tops 
follow parabolic paths. Recent observations, as well as several new, usually 2D, 
simulations have suggested that dynamic flbrils are driven by magnetoacoustic 
shock waves, which originate in the convection zone and photosphere. It has long 
been recognized that along inclined magnetic flux tubes p-mod es leak sufficient 
energy from the gloabal resonant cavity into the chromosphere (IHansteen et al.l 
bOOdlDe Pontieu et all l2007al: iHeggland et alll2007l) . lHeggland et al.l (|2009l) per- 
formed 2D simulations that included the effects of radiation and heat conduction 
and that involved reconnection events induced by waves propagating upward from 
the photosphere/convection z one and were able to produce a number of jet phe- 
nomena among them spicules. iMartmez-Svkora et al.l ( 20091) performed 3D simu- 
lations and showed that some structures having the properties of mottles and DFs 
can be formed by upwardly propagating chromospheric waves through several driv- 
ing mechanisms: collapsing granules, magnetic energy release in the photosphere 
or low chromosphere, p-modes, etc. They suggested that magnetic energy release 
events might be related to magnetic field line reconnection and thus both mecha- 
nisms, i.e. magnetic reconnection and shock waves could co-exist in the jets. Their 
simulations, however, do not reproduce the observed values for the duration and 
maximum height of the jets. Thus, they found durations of 2 - 3 min (instead of 5 
- 10 min) and maximum heights less than 2 Mm (instead of ~ 10 Mm). There is a 
bias also towards rather vertical structures which cannot account for the heavily 
inclined structures that, as it is known, are responsible for the formation of the 
magnetic canopy. Their simulations, on the other hand, showed that half of the 
events they found were driven by magnetic energy release, most likely related to 
the emergence of new field into a pre-existing ambient field. Such conditions, how- 
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ever, do not prevail in plage regions and thus these simulations are still far from 
what works in the real Sun. From the above it seems more plausible to consider 
that both drivers, i.e. waves and magnetic reconnection should be considered as 
responsible for the formation of both mottles and fibrils. 

It is also unclear what the relationship is between the recently named Type 
I and Type II spicules and the traditional spicules. Traditional spicules show mean 
velocities of order 20 - 30kms~^, lifetimes of order 5 - lOmin, and heights of 5 000 
10 000 km. Type I spicules appear to rise up from the li mb and fall back again 



and show similar dynam ical evolution as active region DFs (IHansteen et al.Ll2006c 

De Pontie u et al. i l2007ah . as well as a subset of quiet-Sun mottles ('Rouppe van der Voort et al.l , 

2007 ). For these structu res magnetoacoustic shocks have also been suggested as 

their driven mechanism ((iHansteen et al.l . l2006l : lDe Pontieu et al l l2007b[ l Heggland et al.L 
20071: iMartmez-Svkora et al.l. |2009[ )~ The so-called Type II spicules, on the other 



hand, appear to exhibit upward motion followed by rapid fading without a down- 
ward moving phase. They have shorter lifetimes (10 - 100 s), high apparent ve- 
locities (50 - 150kms~^), and lower widths (150 and 700km). They undergo a 
swaying motion which has been suggested to b e caused by the upward propaga- 
tion of Alfven waves ( De Pontieu et all l2007cl) . Thus it seems that Type I and 



Type II spicules, as seen by SOT on Hinode, look totally different from what 
are considered as traditional spicules described in the literature. It has to be no- 
ticed, however, that the properties of the newly discovered Type I and Type II 
spicules have been derived from observations in the Can H intensity images ob- 
tained with the broad-band filter onboard the Hinode spacecraft, while most of 
the earlier works were based on filtergrams and/or spectra obtained in the Ha 
line. Furthermore, since spicules are so dynamic and many undergo transverse 
motions during their lifetime, in order to infer a reliable measure of their proper- 
ties if a fixed cut almost perpendicular to the limb is used then one has to rely 
only on spicules that did not move transver sely so that their full lifetime can be 
co vered. In a recent w ork Zhang et al. l (|2012l) re-analyzing the same data sets used 



bv iDe Pontieu et al.l (|2007bl) traced the intensity along individual structures in 



quiet Sun and CH. They claimed that they could not find a single convincing ex- 
ample of Type II spicules. Furthermore, they claimed that more than 60% of the 
identified spicules in each region showed a complete cycle, i.e., the majority are 
Type I spicules. Due to these discrepancies, it cannot be clear how the traditional 
spicules are related to the Type I and Type II SOT/Hinode spicules. Further- 
more, what is measured in Can H images is apparent motions and it is unclear 
whether they are associated with bulk mass motions derived sometimes in the Ha 
line. The differences in the heights, lifetimes, velocities and driving mechanisms 
that pointed that there are two types of spicules should also be questioned. Are 
there really two different distributions of the different paramet e rs or one broad 
distribution with a wide range of values? iBrav and LoughheadI ( 1974) gave his- 



tograms of upward and downward velocities in traditional spicules in the range 
±60 - 70kms~^ , lifetimes in the range 2 - lOmin, heights in the range 6-16 Mm. 
iHeristchi and MouradianI ( 1992|) give also a distribution of the apparent velocities 



in the so-called traditional spicules from which it can be seen that they are in 
the range - 70kms~"'^ (see Fig|6l left). Regarding the formation mechanisms, 
[Shibata et al. (1982a) found that: if the explosion (reconnection) occurs below the 
middle chromosphere, the jets are formed as a result of shocks; if the explosion 
(reconnection) occurs above the middle chromosphere the jets are directly driven 
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by reconnection. [Shibata and Suematsul ( 1982b[l using the shock model explained 
why spicules are taher in coronal holes and shorter in active regions. Thus in our 
opinion the question whether there really exist two types of spicules or one type of 
spicules with very broad distributions of heights, lifetimes and velocities, remains 
open. 

Recently an important effort has been undertaken to establish a relationship 
between Type I or Type II spicules and structures on the disk, such as straws and 
RBEs. It is suggested that straws seen in DOT Can H near-limb images appearing 
bright against the very dark background and fast- waving might be identical to the 
fast-waving Type II spicules. Langan gen et al.i (j2008cl) . suggested that Can 8542 A 
RBEs might be the on-disk counterparts of type II spicules observed at the limb 
by the Hinode Can H broad filter. This suggestion was based on the similarity 
in evolution, lifetime, spatial extent and location. RBEs, however, show low blue- 
shifts of 15 - 20kms~"'^, as opposed to the high apparent velocities reported for 
Type II spicules. In order to explain this discrepancy they used Monte Carlo sim- 
ulations and showed that the low observed blue-shifts can be explained b y a wide 
range of spicule orientations combined with a lack of opacity (Langan gen et al.L 
l2008ch . Of course, type II spicules have very short lifetimes, small spatial dimen- 
sions and change very fast a nd thus they are on the limit of wh at can be resolved 
with present-day telescopes. IRouppc van der Voort et al.l ( 20091) from spectral ob- 
servations in the Ha and Can H 8542 A lines reported also that the disk coun- 
terparts of Type II spicules might be the RBEs based on the similarities on their 
appearan ce, lifetimes, longitudina l and transverse velocities and occurrence rate. 
Recentlv. De Pontieu et all ( 201l[) used on-disk Ha observations of RBEs obtained 
with SOT/Hinode and reported that they found direct evidence of a strong corre- 
lation between the RBEs and short-lived brightenings in a wide range of transition 
region (TR) and coronal passbands observed with the Atmospheric Imaging As- 
sembly (AIA) on-board the Solar Dynamics Observatory (SDO) space mission. 
They, also, found evidence that chromospheric spicules observed in Can H with 
SOT/Hinode at the limb in coronal holes are intimately linked to the formation 
of features at (TR) and coronal temperatures and concluded that on-disk RBEs 
and Can H limb spicules are connected. 

It should also be mentioned that in a recent work, lJudge et al.l (|2011 ^) ques- 
tioned whether spicules (and fibrils) are straw-like fine structures (as assumed 
by most authors) or whether they are best described as warps in a sheet. Read- 
ron (2012, private comm.) suggests that some features may appear within «2 
sec, clearly questioning the plasma ejection flux-tube model as this implies Mach 
speeds well in excess of 100. Higher cadence data from present-day imaging spec- 
trometers will allow us to access w hether (or what fra ction) of these features are 
jet/tube features or as suggested bv I Judge et al.l ( 201l|) as current sheets. Spicules 
produced by these current sheets could also be a natural explanation for the Ca n 
H Type II spicules which are very fast and seem to disappear rather than falling 
back to the lower atmosphere. 

Furthermore, it is not known if spicule-like structures seen in UV and EUV are 
some aspects of traditional spicules or totally independent features. Several authors 
suggest that UV spicules may be the result of heating of chromospheric spicules to 
UV temperatures. This would provide an explanation for the reports that spicules 
often "fade from view" when being observed in chromospheric spectral lines (e.g.. 
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iBeckerd (|l968[) . Another possibility is that UV spicules are a transition region 
"sheath" that surrounds the chromospheric spicules. 



10 Summary 

In this review we present observations and physical parameters of fine-scale struc- 
tures observed in the solar chromosphere. The solar chromosphere is dominated 
by a variety of thin, jet-like features omnipresent in quiet, as well as in active 
regions, and on the disk, as well as at the limb. Because of their small size and the 
limitations of ground-based observations, for many decades it was difficult to char- 
acterize completely their properties, and much less to determine with confidence 
their interrelationship and the mechanism(s) responsible for their generation. In 
recent years because of the improved quality of observations acquired by ground 
and space-based instrumentation and, especially, with the current ability to ob- 
serve the chromosphere with high resolution and high time cadence exciting jet-like 
features have been detected that were not recognized earlier. Their properties are 
not yet fully unraveled, while new names have been attributed to distinguish them 
from those previously observed. The nomenclature is nowadays rather confusing, 
but there are strong indications that, at least some if not all of the several observed 
fine-scale structures, are physically closely related, if not the same. 

It seems also that all these small-scale structures are related to the fine-scale 
structure of the magnetic field and its evolution. Magnetic field lines are the most 
likely channels for transporting convective energy to the solar upper atmospheric 
layers. Detailed spatial and temporal high-resolution observations of the magnetic 
fields are crucial for understanding the small-scale structures and investigating 
their role in powering the Sun's outer atmosphere. Magnetic field extrapolations 
are also very important tools in reconstructing the 3-D topology of the upper solar 
layers. Different methods have been proposed so far based on different assumptions. 
Of particular interest is the reconstruction of the 3-D magnetic topology in inclined 
fine-scale structures, such as e.g. spicules. These structures, being magnetic fiux 
tubes following the local inclination of the magnetic field lines, define the layer 
of the magnetic canopy, i.e. the layer where the plasma j3 (ratio of the gas to 
the magnetic pressure) is equal to 1. This layer is very important, since acoustic 
waves generated by the convective motions upon meeting this layer undergo mode 
conversion, refraction, refiection and transmission with important effects to the 
upper solar atmosphere. 

The state-of-the-art in the understanding of small-scale events and their inter- 
relationship requires simultaneous, co-spatial, multi-wavelength, high spatial and 
temporal resolution ground-based observations combined with observations from 
forefront observing space facilities of the same solar region. Such observations will 
permit their follow-up from the photospheric to the lower coronal level of the solar 
atmosphere. They will also permit a better correlation between manifestations of 
the same phenomenon in different atmospheric heights and, hence, are a powerful 
tool for understanding the physics that governs the generation and dynamics of 
fine-scale structures. 

Apart from different types of studies based on analyses of observations, there 
is an on-going development and refinement of non-LTE and cloud model inversion 
codes of spectral line profiles based on better understood physics of spectral line 
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formation and radiation transfer. These codes together with advances in theory 
are an indispensable tool for a more accurate determination of several physical pa- 
rameters that describe the observed fine structures. Quantitative determinations 
of such parameters is an essential ingredient not only for their modelling, but also 
for the modelling of the solar atmosphere as a whole. By their accurate determi- 
nation it will also be possible to define the basic physics of, at least, some of the 
coronal heating processes and the origin and generation of the solar wind. 

These developments go hand in hand with improvements in numerical simula- 
tions and improvement of observational technology. Numerical simulations, using 
3D MHD codes together with forward modelling that produce precise synthetic 
data which are able to mimic accurately improved observations are important to 
correctly interpret them. It is hoped that, in the near future, the combination of 
state-of-the art multi-wavelength observations, along with improved NLTE codes, 
realistic MHD simulations and forward modelling, will provide a breakthrough in 
our understanding of the fine-scale chromospheric structures. 
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